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^ . Abstract: We present a mini-review of the problem of evaluating the energy loss of a 

^ . ultrarelativistic charged particle crossing a thermally equilibrated high temperature e'^e~ 

or quark-gluon plasma. The average energy loss AE depends on the particle energy E 
■ and mass M, the plasma temperature T, the QED (QCD) coupling constant a (a^), and 

the distance L the particle travels in the medium. Two main mechanisms contribute 



> 

, to the energy loss: elastic collisions and bremsstrahlung. For each contribution, we use 

, simple physical arguments to obtain the functional dependence AE(E, M,T,a(^g^, L) in 

different regions of the parameters. The suppression of bremsstrahlung due to the Landau- 
Pomeranchuk-Migdal effect is relevant in some regions. In addition, radiation by heavy 
, particles is often suppressed for kinematical reasons. Still, when the travel distance L 

is not too small, and for large enough energies {E > AP/{aT) in the Abelian case and 
^ , E ^ Mj ^foTg in the non- Abelian case), radiative losses dominate over collisional ones. We 

^ \ rederive the known results and make some new observations. We emphasize, in particular, 

that for light particles {jr? <C ctT"^), the difference in the behavior of AE{E,m,T,a(^g-j, L) 
in QED and QCD is mostly due to the different way the problem is usually posed in 
these two cases. In QED, it is natural to study the energy losses of an electron coming 
from infinity. In QCD, the quantity of physical interest is the medium-induced energy 
loss of a parton produced within the medium. Somewhat unexpectedly, considering the 
case of an electron produced within a QED plasma, we find the same behavior as for the 
medium- induced radiative energy loss in QCD (in particular AE^rad oc at small L), 
despite drastically different behaviors of the photon and gluon radiation spectra, the latter 
being due to the fact that the bremsstrahlung cones for soft gluons are broader than for 
soft photons. We also show that the average radiative loss of an "asymptotic light parton" 
crossing a QCD plasma is similar to that of an asymptotic electron crossing a QED plasma. 
For the radiative loss of heavy particles (M^ ^> aT^), the difference between QED and 
QCD is more evident, even when the same physical situation is considered. 
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1. Introduction 



Evaluating the energy loss of a charged particle passing through usual matter is a 
standard problem in nuclear physics that is very well studied both theoretically and 
experimentally [1]. It is known, for example, that for heavy particles (protons) of 
not too high energy, the main contribution to the energy loss is due to collisions with 
individual atomic electrons, while for light particles (electrons) of similar energies it 
is due to bremsstrahlung. 

The energy loss problem can also be posed for a particle passing a hot ultrarela- 
tivistic plasma.^ In particular, one may ask what happens if a particle carrying color 
charge passes a hot QCD medium. In the limit where the medium temperature T is 
very high, T ^ Aq^^, this medium is a quark-gluon plasma (QGP), i.e., a system 
of quarks and gluons with a small effective Coulomb- like interaction, as{T) ^ 1.^ 
Of course, colored particles do not exist as asymptotic states, but one can imagine 
several thought experiments where the energy lost by an energetic parton crossing a 
QGP could in principle be measured. Consider a fast heavy meson (say, a 5- meson), 
consisting of a heavy quark and a light antiquark, coming from infinity (z.e., created 
in the remote past) and entering a tiny thermos bottle filled with QGP on the lab- 
oratory table, as depicted in Fig. |l|. In the hot environment, the heavy quark sheds 
away its light partner and travels through the plasma losing energy. When it leaves 
the bottle, it picks up a light antiquark or two light quarks to form a heavy colorless 
hadron. The difference in energy between the incoming meson and the outgoing 
heavy hadron roughly coincides with the heavy quark energy loss. 

This situation is probably the cleanest one from the theoretical viewpoint. If the 
heavy quark mass is very large, M 3> T, the density of such quarks in the plasma 
is exponentially suppressed, and thus the passing heavy quark is tagged. But it is 
not possible, of course, to do such an experiment in reality. Hot QCD matter is 
produced in heavy-ion collisions for a few 10^^'^ s, and studying its interaction with 
a i?-meson beam is impossible. It is possible, however, to access the energy loss of 
heavy or light particles produced within the plasma (via a hard partonic subprocess). 
A quark-antiquark pair created with large relative transverse momentum in a hard 
partonic process gives rise to two distinct hadron jets. In heavy-ion collisions, the 
energy loss of the quark between its production and its escape from the QGP softens 
the pr-spectrum of the leading hadrons in the associated jet, when compared to 
proton-proton collisions. This effect, called jet-quenching and first anticipated by 

^In all our study we will consider thermally equilibrated and non-expanding (static) plasmas. 

^We will use the term QGP in this restricted sense which is a natural generalization of the 
conventional definition of a plasma [2]. When E. Shuryak first proposed this name, he thought 
about this analogy [3]. Unfortunately, in realistic heavy-ion collisions, the temperature reached is 
not high enough to have as{T) <^ 1, and whether the system can be reliably described perturbatively 
is questionable. In this paper we will only consider a situation where the effective coupling is small 
and perturbation theory applies. 
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AE = E - E' 



Figure 1: A heavy quark passes a thermos bottle filled with QGP. 

Bjorken [4], has now been well-established by the RHIC experiments [5-8].^ A similar 
effect in cold nuclear matter and at lower energies, namely the attenuation of hadron 
energy distributions in deep inelastic scattering off nuclei, has been observed by the 
HERMES [10] and CLAS [11] collaborations. 

Considering a parton initially produced in a (perturbative) QGP, there are two 
distinct cases. Either the produced parton is light, or it is a heavy quark. For light 
partons of sufficiently high energy, E ^ T, the dominant energy loss mechanism 
(at not too small L) is gluon bremsstrahlung. This problem was earlier considered 
in Refs. [12-22]. To evaluate correctly radiative losses, it is important to take into 
account the generahzation of the Landau-Pomeranchuk-Migdal (LPM) effect [23-25] 
to QCD. In brief, bremsstrahlung is a process where a radiation field accompanying a 
charged particle is shed away. But a newborn "undressed" particle must first "dress" 
with its proper field coat before it can radiate again. If the time needed for such 
dressing {formation time) is large, the radiation intensity and thereby the radiative 
energy loss arc suppressed. For heavy quarks, bremsstrahlung is further suppressed 
compared to the light parton case [26] and the relative contribution of collisional 
losses increases. When the mass of the particle is not too large, M <^ y/aET in 
QED and M <S ^/a^E in QCD, and the travel distance L is not too small, radiative 
losses still dominate over coUisional ones. 

Our goal is to rederive and explain the results in a relatively simple way. We 
will not attempt to perform precise calculations (as far as radiative losses arc con- 
cerned, they are very difficult, maybe impossible to do in a model-independent way) 
and will only give the physical reasoning elucidating the parametric dependence of 

^See Ref. [9] for a recent review on jet-quenching. 
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AE{E, M,T,a(^s), L) in different regions of tfie parameters. This pliysical empha- 
size and an accurate analysis of the radiation spectra for hght and heavy particles 
are some distinguishing features of our review compared to the review presented in 



As emphasized above, studying the energy loss of a parton produced in a QGP is 
of phenomenological interest to heavy-ion collisions. It is, however, also instructive to 
discuss the problem of energy loss in QED. In this case, it is more natural to consider 
an asymptotic (on-shell) particle entering and then leaving a domain containing a 
ultrarelativistic e^e~ plasma, but the production of a QED particle within a QED 
plasma can in principle also be considered. Thus, we found it useful to calculate 
the energy loss of a ultrarelativistic charged particle in all different situations we can 
think of, even though some of them are academic. 

We will first discuss in section ^ the collisional contribution to the energy loss, 
considering the cases of QED, QCD, and of heavy and light particles. For collisional 
losses, the way the particle is produced (in the remote past outside the plasma, or 
initially inside the plasma) is not important. In sections ^ and ^ we discuss the 
radiative energy loss of an asymptotic particle crossing a high temperature plasma. 
Section]^ is devoted to QED, where this physical situation is more natural. We study 
in section ^ the similar problem in QCD, where the "asymptotic" partons should be 
understood as constituents of colorless hadrons when entering the plasma. In section 
H we consider the case of a particle produced in a plasma, both in QED and QCD. 
We show there that the quadratic dependence of the induced radiative loss on the 
plasma size L at small enough L is not a feature specific to QCD. It also holds in QED 
with the same parametric dependence (up to logarithms). The differential gluon and 
photon radiation spectra are qualitatively different, however. Since it might be useful 
to jet-quenching phenomenology, the induced gluon energy spectra corresponding to 
light and heavy quark radiation are discussed in some detail. Finally, we briefly 
summarize and give some general remarks in section ^. 

2. Collisional energy loss of a fast charged particle 

The first calculation of the collisional loss of a fast charge crossing a hot plasma 
is due to Bjorken [4]. This was done in the context of QCD, and was the basis of 
Bjorken's proposal to use jet-quenching as a signature of the QGP. Bjorken's result 
for the collisional loss (per unit distance) of an energetic light parton (light quark q 
or gluon g) reads 



■^Think for instance of the energy loss of a charged lepton produced in a heavy-ion coUision 
(though it would be small in this case, due to the smallness of the QGP size compared to the lepton 
mean free path). 



Ref. [27]. 
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Figure 2: The typical graphs for collisions of a muon with plasma particles. 

where Uf is the number of thermally equilibrated quark flavours, Cr = Cp = 4/3 
(quark) or Cr = Nc = 3 (gluon), and /i is an effective infrared cut-off. To the 
logarithmic accuracy, to which is derived, /i can be taken as the Debye screening 
mass in the QGP, fi ~ gT. 

More detailed studies of the collisional losses in QED and QCD plasmas were 
performed for instance in Refs. [28-34] . We review below this problem, emphasizing 
the difference between the two cases of a tagged (heavy) or untagged (light) particle. 

2.1 Hot QED plasma 

2.1.1 Ultrarelativistic muon 

We consider the case of a ultrarelativistic muon of mass M and four-momentum 
P = (EjP) passing a hot e~^e~ plasma of temperature T, with E ^ M ^ T. The 
second inequality ensures that there are no muons in the heat bath. The muon can 
lose its energy in either Coulomb collisions with electrons and positrons (Fig. |^a), or 
Compton collisions with photons (Fig. 0d). 

Consider first the losses due to Coulomb scattering. The differential Coulomb^ 
cross section is of the form 

dcTCoulomh ^2 2^ 



2^2 



dt {t — 1^ 

where t = = {K — K'Y is the Mandelstam invariant momentum transfer, and /i ~ 
eT is the Debye screening mass in the QED plasma. The total Coulomb scattering 
cross section reads 

(^Coulomb ~ ^ ~ ™ . (2.3) 

The momenta K and K' in Fig. |^ refer to thermal particles and are of order T. 
Note that they are slightly off mass shell, ~ K'^ ~ /i^, due to medium effects, 
which does not affect, however, the estimates below. Since the integral for the total 
collisional loss is saturated by the region |t| ^ /i^ (as we will see in a moment), we can 



^Following the usage adopted in the Hterature on this subject, we use the word "Coulomb" in a 
generahzed sense, inchiding also screened Coulomb, i.e., Yukawa. 
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write t = -2KQ = -2{\K\Qo- K -Q). We have in average Qo t/{2\K\) ~ -t/T. 

The mean energy loss in a single scattering is thus 

(AE)iscat. ^— [ (it^^^^^~aTln%i. (2.4) 

^Coulomb J at I fl 

The logarithm arises from the broad logarithmic interval /i^ ^ \t\ <^ |t|max, imply- 
ing (to logarithmic accuracy) that the energy transfer in a single collision is small 
compared to E. 

Introducing the Coulomb mean free path, 

_ 1 1 

^Coulomb = ~ 5 (2.5) 

'T-CCoulomb Oil 

where n ~ is the particle density in the plasma, we estimate the rate of energy 
loss per unit distance as 

1 scat. 2 ATI 2 1 '' max /p, «\ 

~ — ~ a T m — — . (2.6) 

ax Acoulomb 

The maximal transfer |t|max is given by 

I.I - Mr E^T^ 

|t|max- ^ ^ M^ + 0{ET)' ^ ^ 

where we used s = {P + K f = + 2P ■ K = + O (ET). Thus, we see that 
the logarithmic factor in ( p.6|) takes a different form in the mass domains -C ET 
and M2 > ET, 



In ^ (M^ < ET) 



'^-^Coulomb 2t-i2 



dx 



a'T' { . (2.8) 

21ng (M2»i?T) 



Two remarks are in order here. 



(i) Strictly speaking, using (|2.2|) for the differential cross section is not correct. 
First, it is written for static and scalar plasma particles, while the particles 
are moving and have spin. Second, the graph in Fig. ^ was evaluated with 
the model expression l/{t — fi^) for the photon propagator, while the actual 
expression is more involved. Third, neither S'-matrix nor cross section is well 
defined in medium, because there are no asymptotic states, and the relevant 
quantity is not 0"*°*, but the damping rate C of the ultrarelativistic collective 
excitation with muon quantum numbers, related to the muon mean free path 
by A = C^^. This damping rate was evaluated in Ref. [35], with the result^ 

-1 



^aTln (C/a)) 



[2.9] 



^The calculation was done for QCD, but the Abelian result is directly obtained from the result 
for (^q by setting Ci? = 1 in Ref. [35]. 
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Thus, the mean free path involves an extra logarithm of the coupling constant 
as compared to our estimate (^^ ). This is because (|2.9| ) depends on mag- 
netic interactions in addition to Coulomb scattering (see also Ref. [36]). Thus, 
strictly speaking, it is ( |2.9| ) which should be taken as the definition of A to 
be used in sections ^ to ^. However, keeping track of the logarithms of the 
coupling constant is a difficult problem, which we will not address. The only 
logarithms we will keep are those depending on the particle energy E. Hence, 
in all our study we will use A ~ l/(aT), which coincides with the Coulomb 
mean free path 



(ii) A more standard way to define the mean free path is not as in ( |2.5[ ), but rather 

by A*'' = 1/(?7.(T*''), where 



a'' = J da {1- cose) (2.10) 

is the transport cross section involving an additional suppression factor for 
small angle scattering. The transport mean free path A*"^ conveniently de- 
scribes standard transport phenomena associated with collisions (viscosity, elec- 
tric conductivity, ...).'' For the problem of collisional energy loss, the use of 
the definition ( p.lO| ) is equally warranted, when deriving ( |2.6| ) (the estimates 
for {AE)i scat and A would be different, but their ratio would not change). 
The scale (^3]) appears to be rather elusive also for radiative losses. Almost 
all results there depend not on A as such, but on the combination q = /i^/A, 
which is a transport coefficient. The "isolated" scale A only enters arguments 
of certain logarithms. We will come back to the discussion of this point, when 
concluding in section ^. We only stress now that, irrespectively of whether it 
is observable or not, the notion of mean free path as defined in (|2.5| ) or ( |2.9| ) 
proves to be very convenient and instructive, and we will stick to this definition 
throughout the paper. 

Let us now turn to collisional losses due to Compton scattering (Fig. ^d), and 
focus on the region -C ET. The differential Compton scattering cross section 



reads, for < |m| < s ~ ET, 



dt su 
This gives the total Compton cross section 



C^O'Compton _ (2 11) 



s ET 

^Compton ~ V M2 ^ ET JP' ^^'^^^ 



^For a thorough discussion of transport phenomena in a QGP, see [37] 
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arising to logarithmic accuracy from the domain <^ |m| <^ s ~ |t| (recall that 
s + t + u = 2M^) 
scattering, 



Similarly to (|2.4|) , we find for the energy loss in a single Compton 



1 



1 Compton scat. 



7. '^'^Compton ^ 7-T 
at ^ _ ^ ^ . 

dt T T 



(2.13) 



^Compton 

Thus the characteristic energy transfer is of order E. (The experimental projects to 
produce energetic photons by scattering energetic electrons off laser beams are based 
on this property of Compton scattering [38]). Introducing the Compton mean free 
path 



A, 



1 



E 



Compton 



^'^Compton 



> A 



Coulomb 



we obtain 



dE< 



Compton 

dx 



E 



ET 
M2 



(2.14) 



(2.15) 



Acompton 

Note that the logarithm in ( p.l5| ) arises from the same logarithmic integral (over u) 
as for the total Compton cross section ( |2.12| ), and that it is present only in the mass 
domain < ET. 

We see that the losses due to Coulomb and Compton scattering are of the same 
order, although the two processes differ drastically. Compton scattering is rare, 
Acompton ^ Acouiomb, but, as mentioned above, it is very efficient in transferring 
energy. Summing the Coulomb and Compton contributions to the collisional loss in 
the domain < ET, we get^ 

dEcow 



dx 



TC 



ET lET , , 

1.^ + 51.5^ + 0(1) 



(M^ < ET) . (2.16) 



For Compton scattering, when <^ ET we have \u\ <^ s to leading logarithmic 
accuracy, implying that \t\ ~ |t|max — s. Thus, in this mass region Compton scat- 
tering is characterized by a final state consisting of a soft muon and a hard photon 
(of energy ~ E) ejected from the plasma. This is to be confronted with the charac- 
teristic Coulomb scattering process, where the energy transfer is small. Beyond the 
leading logarithm, configurations where the final muon and the scattered thermal 
particle share similar fractions of the initial energy E also contribute to dE/dx, for 
both Coulomb and Compton contributions [33]. 

Finally, let us note that for ^ ET, the Compton logarithm in ( p. 161) should 
be dropped, and the Coulomb logarithm is modified, see ( p.8|) , giving 



dE. 



coll 



dx 



27r 



a'T' 



T^T 



(M^ > ET) . 



(2.17) 



^See Ref. [33] for the exact calculation, where the constant beyond logarithmic accuracy is also 
evaluated. 
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2.1.2 Ultrarelativistic electron 

It is tempting to estimate the coUisional loss of an energetic electron crossing an 
e~^e~ plasma by replacing in ( |2.16| ) the muon mass M ^ ^ by the electron thermal 
mass in the medium mth ^ eT fi. However, as noted above, the leading Compton 
contribution corresponds to the situation where the incoming particle loses almost 
all its energy. Thus, it seems impossible to distinguish the final soft electron from 
a thermal electron. In other words, when an energetic electron becomes soft after 
interacting with the plasma, it effectively disappears. In addition, the incoming 
electron can be annihilated with a thermal positron. A similar situation arises when 
discussing positron energy loss in usual matter. 

One way to better define an observable energy loss in this case is to require the 
final electron to be hard enough, say, with energy E' > E/2. This constraint makes it 
very unlikely that the final energetic electron is a thermal electron, due to the expo- 
nential suppression of the Fermi-Dirac thermal weight, np{E') ~ exp {—E'/T) -C 1 
for E' > E/2 ^ T. Demanding the presence of an energetic electron in the final 
state allows one to discard the annihilation channel, as well as the leading (logarith- 
mic) Compton contribution. Only the Coulomb contribution, which to logarithmic 
accuracy corresponds to small energy transfers, should be kept. 

With this setup, we obtain from ( |2.16|) 



dE, 



coll 



dx 



In— + 0{1] 



(2.18) 



-,E'>E/2 '-' 

2.2 Quark gluon plasma 
2.2.1 Tagged heavy quark 

Here we consider the case of a ultrarelativistic heavy quark {E ^ M ^ T) crossing 
a QGP. For clarity we focus on the limit < ET. (As in QED, for > ET the 
Compton logarithm should be dropped in what follows, and the Coulomb logarithm 
modified, see ( p.8|) .) Purely coUisional energy loss does not depend much on the 
production mechanism of the heavy quark. The latter can be produced in a hard 
process within the medium (like in heavy-ion collisions) or preexist within a heavy 
meson coming from infinity. The crucial difference between these two situations 
will reveal itself when studying the radiative energy loss of a quark induced by its 
rescattering in the plasma. 

The calculation of the heavy quark coUisional loss in a QGP is similar to the 
case of a muon crossing an e~^e~ plasma. The main change consists in the running 
of the coupling a^. This is essential to take into account, not only to improve the 
accuracy of predictions, but also to obtain the correct energy dependence of dE /dx. 

The contribution to dE/dx from Coulomb scattering of the fast heavy quark 
off thermal quarks and gluons is easily inferred from the QED case. The Coulomb 
differential cross section is oc a^, and the scale at which to evaluate is given by 
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the invariant momentum transfer t itself. Due to the running of a^, the logarithmic 

integral appearing in the fixed coupling (QED) expression (|2.4|) is thus modified 
to [32] 

, f^^ d\t\ d\t\ 2, , , , 

« / -Uo^lit)- (2.19) 

Using aait) ~ 1/ In (|t|/A^), the r.h.s. of the latter equation can be exactly inte- 
grated, and we can rewrite it as 



a In 



ET 
1^ 



as{fi'^)as{ET) In 



ET 
1^ 



(2.20) 



A similar discussion applies to the contribution from Compton scattering off 
thermal gluons.^ To logarithmic accuracy, Compton scattering is dominated by u- 
channel exchange. The relevant scale determining the coupling in the differential 
cross section for this contribution is ~ 0{u). The total Compton scattering cross 
section in QCD is obtained from the QED expression (|2.12|) by replacing 



a 



In other words, 



a In ■ 



ET 



Af2 



ET 



d\u\ 



\u\ 



ET 



M2 



d\u\ 



\u\ 



al{u) . 



as{M'^)as{ET) In 



ET 



(2.21) 



(2.22) 



Using (|2.20|) and (|2.22|) in ( p. 161) , we obtain for the fast heavy quark collisional 
loss in the limit <C ET (after performing the thermal average over the target 
quarks and gluons and introducing color factors [34]), 



dE, 



coll 



dx 



(l + ^) «,(/x2)a,(ET)ln 



ET 
1^ 



+ ^a,iM')as{ET)\n^ + 0{al) 



(2.23) 



The term beyond logarithmic accuracy ~ O (a^) was determined in Ref. [34]. Simi- 
larly to the QED case, the Compton leading logarithm in (|2.23|) corresponds to final 
state configurations with a soft (but tagged) heavy quark jet and a hard jet initiated 
by a gluon of energy ~ E knocked out of the plasma. 



^In QCD, the terms Coulomb and Compton refer not to different processes, as tliey do in QED, 
but to different kinematical regions associated with the same process. Thus, the amphtude M (Qg 
Qg) is dominated by the Coulomb diagram with soft gluon exchange when \t\ is small. The same 
amplitude is dominated by the Compton-like diagram corresponding to u-channel exchange when 
\u\ is small. 

-^°We note that such configurations are presently not counted in the RHIC experimental setup, 
due to a lower energy cut-off used in the selection of heavy quark tagged events. Under those 



experimental conditions, only the Coulomb leading logarithm should be kept in (2.23) 
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2.2.2 (Untagged) light parton 

Similarly to the case of QED, it would be misleading to pretend obtaining the energy 
loss of an energetic light parton by replacing in ( p.23| ) the mass M by the parton 
thermal mass ~ gT. First, even though this is not done in practice, it is theoretically 
possible to observe the events with a soft tagged heavy quark. But for a light (and 
hence untagged) parton, this is impossible. Second, light quarks may annihilate with 
the light antiquarks in the heat bath, and this adds to the intrinsic uncertainty of 
what energy loss of a light parton is. 

As far as Compton scattering is concerned, the situation is even worse than 
in QED, where the detection of an energetic photon in the final state would at 
least signal that a Compton scattering occurred. In QCD, it is very difficult to 
distinguish the hadron jet initiated by a hard final gluon produced in Compton 
scattering from the hadron jet initiated by a hard final quark having undergone soft 
Coulomb exchanges. 

For a light (or, more generally, untagged) parton, the observable energy loss must 
be defined, at the partonic level, with respect to the leading {i.e., most energetic) 
parton. When E' < E/2, or equivalently \u\ < s/2, the corresponding energy loss 
is thus AE = E - \K'\ = E' - \K\ ^ E' ^ \u\/{2\K\). When E' > E/2, we 
have AE = E — E' ^ |t|/(2|i^|), as in the case of a tagged parton. The Compton 
contribution to dE/dx is thus of the form 



dE, 



Compton 

dx 



n 



d\u\ 



M2 



a: 



s\u\ 



'Ui\ \f\ 
^-^Q{s/2-\u\) + ^Q{\u\-s/2) 



(2.24) 



We observe that, in contrast to the tagged case, the Compton contribution has 
no logarithmic enhancement. The reason is clear: the logarithm ~ In (E'T/M^) in 
( p. 23 ) came from the small u region. But with the "untagged" physical definition 
( p.24|) , small u corresponds to small energy transfer, and does not give any important 
(logarithmic) contribution to the energy loss. Within such a definition, we also note 
that the annihilation channel contribution does not yield any logarithm either. 

Thus, as far as leading logarithms are concerned, we are left with only the t- 
channel Coulomb logarithm arising from the broad interval /i^ -C |t| <^ ET . This 
Coulomb logarithm can be read off from ( |2.23 ). Generalization to the gluon case is 
obvious and we obtain 



dE, 



coll 



dx 



CrtcT' 



q.g 



1 + 



ET 



as{fi^)as{ET) In— + 0{al) 



(2.25) 



with Cr = 4/3 (3) being the quark (gluon) color charge. Replacing as{fi'^)as{ET) 



a 



2 ;v 



in ( |2.25| ), we recover the fixed as Bjorken's result (|2.1| ). We stress that Bjorken's 



result is valid in the logarithmic approximation and specific to the untagged experi- 
mental setup. 
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3. Radiative loss of a fast asymptotic charged particle — QED 



We will now discuss the radiative energy loss due to bremsstralilung. In this section, 
we consider the case of an asymptotic charged particle {i.e., produced in the remote 
past) crossing a plasma layer of finite size. This is a natural experimental setup 
in QED, where it is possible to prepare an on-shell energetic electron entering the 
plasma with its already formed proper field "coat" . We thus consider the QED case 
first. In section ^ we will study the somewhat academic but instructive case of an 
"asymptotic quark" crossing a QGP. 

3.1 Fast electron 

Here we evaluate the radiative energy loss of an on-shell energetic electron going 
through a ultrarelativistic e~^e~ plasma. We will assume E ^ T and fi ^ eT ^ m, 
where m = rUf. is the electron mass. The electron is scattered by the plasma particles, 
changes the direction of its motion and emits bremsstrahlung photons. 

L <^ X: Bethe-Heitler regime 

Let us first consider the case of a very thin plasma layer of size L -C A, where the 
electron mean free path A is the characteristic distance between subsequent elastic 
scatterings.^^ As discussed in section 0, A is given by (|2.9|) or rather by (|2.5| ) since 
we neglect logarithms of the coupling constant in our study. The probability that 
the incoming electron undergoes a Coulomb scattering is ~ L/\ ^ 1, and the 
probability to have several scatterings is further suppressed. The average energy loss 
after crossing the length L is thus 

AE(L « A) ~ ^AE-d,, , (3.1) 

where /S.Ell^^^ is the radiative energy loss induced by a single Coulomb scattering. 
It is obtained from the photon radiation spectrum derived by calculating the two 

^^In our paper, L denotes the distance travelled by the particle in the plasma, to be distinguished 
from the plasma size Lp. For thermal equilibration to occur, the size of the medium should be 
much larger than the transport mean free path of plasma particles (which is of the same order 
as the transport mean free path of an energetic electron), Lp ^ A*'' ^ l/(a^T) ^ A. Thus, the 
situation i <C A, implying L ^ Lp, is quite unrealistic for an asymptotic electron, for which we 
expect L ^ Lp, except if the electron is crossing the plasma near its edge. (It is somewhat more 
natural to consider the limit L <C A in the case of a particle produced within a plasma, as we will 
do in section ^.) 

However, in order to better understand what happens in the more physical situation with larger 
L, we find it instructive to consider first the case L <C A. 

^^The radiative energy losses induced by a single scattering will be referred to as Bethe-Heitler 
(BH) losses in our paper. 
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Figure 3: The two diagrams for plioton radiation induced by electron Coulomb scattering. 



diagrams of Fig. ^ 



rad 
scat. 



dL 



rad 



uo duj . 



(3.2) 



The result for the radiation spectrum does not depend much on the nature 
of scatterers. It depends mainly on the characteristic scattering cross section and 
thus on the characteristic momentum transfer q^. As discussed previously, we may 
consider only Coulomb scatterings. Thus, we can replace the plasma particles by 
static sources and use the form (|2.2|) for the differential cross section. The integral 



giving the total Coulomb cross section ( ^.3[ ) is saturated by the values |t| ~ ~ /i^ 
and, in most estimates (with few exceptions to be spelled out later), we can assume 
to be on the order of the Debye mass /i. 

In the soft photon approximation uo = \k\ which is sufficient for our 

purposes,^" 



the photon radiation intensity reads 



dL 



rad 



i=l,2 



P ■ Si P' ■ Si 



P-k P' -k 



d^k 



(27r)32w 



(3.3) 



where Si are two transverse photon polarization vectors. When the photon radiation 
"angle" = k±/u! and the electron scattering "angle" 9s = q±/E are small, 9, 9s -C 1, 
the sum over photon polarizations gives 



dIr!>A — — d^9 J?' 



J. 



9' 



QI2 Q2 



9^ + 6*2 



9' 



+ 9l){9'^ + 9: 



9r, 



9„ 



02 + ei 9'^ + 91 



(3.4) 
(3.5) 
(3.6) 



-'^^In fact, the characteristic frequency of emitted photons contributing to AE,-ad is of order E 
(see e.g. the spectrum ( |3.9D below). However, this modifies the result obtained in the soft photon 
approximation only by numerical factors, which is not our concern here. 
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where 9m. = m/E and 9' = 9 — 9s. 

We will be interested in the small and large mass limits, corresponding respec- 
tively to 9m <^ 9s and 9^,^ 9s- Thus, we quote: 



/ d?9 -»r, , 9l c^-^rad 2a . 9'^ 

9m^9s^ / — j/^ln^jf ^ c^-^^-ln^, (3.7) 



m 



(3.8, 

y 277 3 6'4 duj 37r 6'4 

Eq. ( p.7|) is valid to logarithmic accuracy. The logarithm arises from the first term of 
(|]6D and from the angular regions 9m<^9 <^9s and 9m<^9' = \9 -9s\ < 6*3. The 
asymptotic expressions ( p.7|) and (3^) are conveniently incorporated in the following 



interpolating formula for the energy spectrum, 

d/,ad 2a f d'9 -.2 2a / 9's \ , ^ ^ ?! A 

The exact expression displays the same qualitative behavior, but is more complicated 
[39]. 

When q± m (implying 9s 9m), we obtain 

/E JT 2 
-^uduj^aEln^^, (3.10) 
dio m'^ 

and for the Bethe-Heitler (BH) energy loss per unit distance, 

^(L«A) ~ -a^i^Tln^ • (3.11) 

dx A 

Let us now discuss a possible modification of this result due to Compton scat- 
tering. Essentially, there is none. (This is in contrast to the coUisional loss. We have 
seen in section 2 that the contributions due to Compton and Coulomb scattering are 
of the same order there.) Indeed, for Compton scattering, ^El^^^ is of the same 
order as for Coulomb scattering, ^E^'^^^ ~ aE, but the corresponding mean free 
path ( |2.14| ) is much larger and 



J TP raid rp 

9^^(L « A) ~ a'T' (3.12) 

^Compton 



is much smaller than the BH contribution (p.ll|) induced by Coulomb scattering. 



(The Compton contribution ( |3.12 ) is even smaller than the coUisional loss ( 2.18| ) 



and will not be mentioned any more in this paper.) 

We see that, when L <^ A, the BH radiative loss ( p.ll| ) of an energetic electron 
crossing a ultrarelativistic plasma is much larger than its coUisional loss (f2.18| ), 

L«X - |»1. (3.13) 
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It is instructive to compare this with the energy losses of a ultrarelativistic par- 
ticle in usual matter [1] , consisting of nonrelativistic electrons and static nuclei. The 
basic difference between usual matter and an e+e" plasma is that, in the former, 
screening effects do not play an important role. For sure, the electric fields of indi- 
vidual electrons and nuclei are screened at atomic distances ~ l/(am), but these are 
comparatively large distances (see the footnote below). Screening in usual matter 
can affect the arguments of certain logarithms, but is otherwise unimportant for or- 
der of magnitude estimates. In the following digression on usual matter, we neglect 
logarithms. 

Consider first coUisional losses. Weighing the e^e^ —>■ e^e^ Coulomb differential 
cross section da/dt ^ /t^ by the energy transfer AE(t), we get 

(7 T-i \ usual matter n 2 

where n is the number of atoms per unit volume and Z is the number of electrons in 
an atom. Using t = -2K ■ Q ~ -2{mQo - K -Q), we have Qq = AE{t) ~ \t\/{2m). 
We obtain, up to some logarithm, 

^^ysual matter ^^^2 

Let us note that the order of magnitude of collisional energy loss in a hot plasma is 
recovered from ( |3.15| ) by replacing Z ^ 1, n ^ and m ^ T. 

Now we discuss radiative losses. Let us assume (as we also did for a hot plasma) 
that the matter layer is thin enough, such that AE{L) <^ E. In addition, let 
us first assume the scattered energetic particle to be an electron and the medium 
hydrogen [Z = 1). Then the radiative losses due to scattering by electrons and nuclei 
(protons) are of the same order. When q± <C m, the characteristic energy loss in a 
single scattering is ~ aE times the suppression factor ~ q\/rn? ~ |t|/m^, see (|3.8| ). 
Hence, instead of ( p. 141) and (|3.15|) we obtain for the BH radiative loss 

hydrogen .2 u| n E , , 

dx J -QYi J ""^ 

Neglecting logarithms, the estimate ( |3.11| ) for radiative losses in a ultrarelativistic 
(thin) plasma is obtained from ( p.l6| ) by replacing n and m —>■ fi eT, instead 

of m — >^ T as we did in the case of collisional losses. 

For usual matter, the characteristic ratio of radiative and collisional losses is 

^^^^ (3.17) 



dEcoii m 



^^This can be understood by noting that in the plasma case, because of screening at the scale 
fi, only the region \t\ ^ fi^ ^ contributes. (In usual matter, screening occurs at some scale 
l^lmin ^ m? ■) Thus, the factor \t\/m?' in ( 3.16| ) disappears and the t-integral is saturated by |t| ~ ^ 



2 
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with an additional suppression factor a compared to (|3.13|) , due to the different 
screening properties of plasma and usual matter. That is why radiative losses in 
usual matter dominate only at comparatively high energies E ^ mja. For electrons 
in hydrogen, the critical energy where the radiative and coUisional losses are equal 
is ~ 350 MeV [40]. 

The physics of coUisional and radiative losses differ in one important respect. The 
energetic particle loses a tiny fraction of its energy in a collision with an individual 
electron, much like a cannon ball loses a tiny fraction of its energy in a collision 
with an individual air molecule. The drag force dp/dt = dE/dx is an adequate 
physical quantity to describe this. On the other hand, about half of the original 
particle energy is lost during a single radiation act. The estimate ( |3.16| ) refers to 
an average drag force, while the fluctuations of this quantity are very large. That is 
why radiative losses are usually not described in terms of the drag force ( |3.16[ ), but 
in terms of the radiation length Xq - the average distance at which about half of the 
energy (more exactly, the fraction 1 — 1/e) is lost. It is especially sensible bearing in 
mind that —dE-Bu/dx oc E and the energy thus decreases exponentially. 

We will see later that, though the radiative energy losses of light partons passing 
a QCD plasma strongly fluctuate by the same reason as electron radiative losses, the 
notion of radiation length is not convenient there, because the linear BH law ( p.l6| ) 
is not realized and the energy dependence of dE/dx is more complicated. We can 
also mention right now that the radiation spectrum of heavy enough quarks (but 
not of heavy Abelian charged particles!) turns out to be soft, so that the physics 
is more similar to the physics of coUisional loss and the drag force fluctuations are 
suppressed. 

When Z > 1 and the incoming particle is heavy, M ^ m, two new effects come 
into play. First, radiation mainly occurs when the incoming particle is scattered on 
a heavy nucleus, due to the enhancement factor (square of the nucleus charge) 
compared to the hydrogen case. The coUisional loss is still due to scattering off 
electrons and is only enhanced by the number of electrons Z. Second, the radiation 
intensity is suppressed at small momentum transfers by the factor ~ q^/M"^. We 
obtain: 

dE^n ZaEm 

dE-^r ^ ■ ^'-''^ 

The suppression factor q^/M"^ ~ fi'^/M'^ is effective also for radiative losses of 
a massive particle passing through a ultrarelativistic QED plasma (see ( |3.8| ) and 
( |3.37| )). The meaning of the condition /i ~ eT ^ m imposed at the beginning of the 
section has now become clear. When fi ^ M, a particle of mass M passing through 
a plasma can be considered as light. When M ^ yU, it can be considered as heavy, 
and its radiative loss in a QED plasma is suppressed (at least for L ^ A, see ( p. 37] ) 
below) by the factor^^ ~ /i^/M^. 
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This suppression, which is due, as can be inferred from (3.8), to the suppression of radiation in 
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L:$> L*: LPM regime 

After this digression we come back to the case of a hot plasma, and consider a 
medium of very large size. The estimate (|3.1| ) is correct for small L <^ A, the factor 
L/X bearing the meaning of the electron scattering probability. One can be tempted 
to extend it also to the region L S> A, with the factor L/X being interpreted as the 
number of electron scatterings, but this is not correct because the basic assumption 
on which it is based, namely, that a photon is emitted in each scattering, does not 
hold. 

Indeed, the radiation process takes time. The particle is not ready to radiate 
again until it grows its accompanying coat of radiation field. Strictly speaking, it 
takes an infinite time for the coat to be fully developed such that the particle can be 
treated as an asymptotic state. But if we are interested only in emission of photons 
in a particular wave-vector range, we are not obliged to wait eternity and should only 
be sure that the corresponding Fourier harmonics of the radiation field are already 
present in the coat. The length at which a harmonic corresponding to radiation of 
the photon of frequency u at angle 6 is formed, the formation length, reads^^ 

(3-19) 

where ~ 2ii^cij(l — cos 6) ~ Eu)6'^ is the virtuality of the internal electron in Fig. 

The formation length ~ 1/{uj6'^) can be interpreted as the length at which a 
photon of energy u emitted at angle 6 acquires a phase of order 1 in the frame 
moving with the particle. The latter condition is, indeed, 

(phase disbalance) ~ cut - k\\i ~ i{u - Vto"^ - ujW) ~ iuO^ ~ 1 , (3.20) 

where we assumed that the electron moves with the speed of light, t = i (one can do 
it as long as k± c:^ uj6 ^ m). 

The crucial observation is that, even if the electron is scattered several times 
before travelling the distance if{uj,6), it cannot emit (to leading order in a) more 
than one photon of energy uj at angle 6. Indeed, in this case one cannot talk about 
independent photon emissions, and we are dealing with coherent emission of a single 
photon in a multiple scattering process. As a result, photon emission with £f{uj,9) 3> 
A is suppressed compared to the situation where it would be additive in the number of 
elastic scatterings. This effect is known as the Landau-Pomeranchuk-Migdal (LPM) 
effect [23-25]. 



the cone 9 < 6m = M/E, is sometimes called dead cone effect [26]. 

"'^^For om pmposcs, it is better to talk about formation length rather than formation time. In the 
rest frame of the virtual electron, the formation time/length is of order In the laboratory 

frame, it is multiplied by the Lorentz factor E/V^. In general, the angle 6 appearing in ( 3.19| ) 



is the angle between the photon and the electron from which it is radiated (the final electron in 
Fig. |a). 
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At fixed uj and when m, the radiation spectrum ( p.9|) induced by a single 
scattering arises from the angular domains 6*^ -C 6* -C 6*^ and dm <^0' = \6 — 6s\ -C 6*^, 
i.e., from formation lengths 

<^if{uj,e)<^^. (3.21) 



In the integrated spectrum, we have u) E, and thus the BH radiative loss ( p.ll| ) 
arises from photon formation lengths 

^<^if{E,9)<^^. (3.22) 

For E ^ T, we have E/fi"^ ^ A, and thus the radiated photons contributing to the 
BH loss (|3.11|) are formed far away from the medium layer of size L <^ X. 

What happens when the medium size increases, L ^ A? Since, for a thin medium 
of size L ^ A, the typical formation length of radiated photons is larger than E/fi"^, 
we would naively expect the energy loss to be approximately given by AEl'^^^^ (see 
( P.10| )) as long as L < E/ fi"^, the entire medium acting as a single effective scattering 



center. However, this is not so. The critical size L* beyond which AE{L) strongly 
differs from AE^'^^^^ happens to be much smaller than E/fi"^. It is approximately 
given by the geometric average of A and E/fi^, 




A < L* ~ J J < (3.23) 



To find how this scale arises, consider the limit L ^ oo. In this case, all radiated 
photons are formed in the medium, and the estimate of the typical photon formation 
length i^'^'^ is modified compared to the vacuum case. To see that, substitute for 6"^ 
in the estimate (^.19]) the typical electron deviation angle 6*^ at the length if'^'^, 



/imed 



(3.24) 



where N = P^'^^ j A is the number of electron scatterings on the distance i™'^^. The last 
estimate in ( p.24|) relies on the Brownian motion picture, namely on the assumption 
that the transverse momentum transfers in successive elastic scatterings are not 
correlated. We obtain 

Substituting ~ i?, we thus get the estimate ( p.23| ) for the characteristic in-medium 
formation length L*. 

Scattering centers located within a distance ~ L* act as one effective scattering 
center, inducing single photon radiation. For L ^ L*, the medium contains L/L* 



r'H-W^- (3-25) 
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effective scattering centers, and the estimate for tfie radiative energy loss is obtained 
by multiplying this number by the energy lost in single photon radiation/'' 



A^(L>L*) ~ aE^ ~ aL^^E ~ L VET'^ . (3.26) 

This estimate has been obtained by assuming that the typical momentum trans- 
fer squared after electron scatterings is q'^{N) ~ N^"^ (/x being the typical transfer 
in one scattering). In fact, this is correct only when the g_L-distribution in single 
elastic scattering decreases rapidly enough at large q±, such that the average (g^) is 
well-defined. However, for Coulomb scattering, the integral 

iA) = I^AAj^r^ (3.27) 

diverges logarithmically when the ultraviolet cut-off fzilj^ax ~ I I max ~ ET oo. 
It is possible to show (see Ref. [14], or Appendix A for an alternative derivation) 
that the typical transfer after A^ Coulomb scatterings scales as ( A^ In A^) /i^ instead 
of Nfi'^. The presence of this logarithm somewhat modifies the estimate (|3.23| ) of L*, 
but otherwise does not affect the above heuristic derivation. We thus find: 



/2 



XE XE , , 

(3.28) 



fiHn{L*/X) \j fiHn{E/{Xfi^)) 
AE(L>L*) ~ aE^ ~ a^L^ET^ \n{E/T) . (3.29) 

The parametric dependence of the latter result agrees with that found in Ref. [41]. 
At both small L <^ A (see (|3.1|) ) and large L ^ L* (see ( |3.29|) ), the radiative energy 



loss is proportional to L. But the proportionality coefficients in those two regions 
are different. At large L, the slope is smaller as a result of LPM suppression. The 
behavior of AE{L) is represented schematically in Fig. ^. 

It is worth noting that for L ^ L*, the in-medium energy spectrum of the 
radiated photons can be easily obtained from ( |3.25| ), 



where we introduced the energy scale 



l2 2 

coc ~ ■ (3-31) 



"'^^This energy is given by (3.10), but without the mass-dependent logarithmic factor. This is 
because our estimate (f'^' ^ ^ L* displays no logarithmic spread, contrary to ( 3.22 ). This can also 



be expressed (using ( 3.19| )) by noting that the electron travelling in the medium is off mass shell, 



with a characteristic virtuality of order nf^^^^ ^ E/L* ~ ^ Efi'^/X ~ aVET^. This is in contrast 
with the logarithmic spread ^ <C /i^ in vacuum. 



- 19 - 




Figure 4: Schematic plot of the radiative energy loss of an asymptotic light (m <C //) and 
fast {E S> T) charge crossing a hot e~^e~ plasma, as a function of the distance L travelled 
in the plasma. The dotted line represents the collisional energy loss AEcoii{L) ~ a 



Including the correct logarithmic factor as discussed above, the spectrum (|3.30D reads 



Integrating the latter spectrum up to ~ E, we recover (|3.29|) . The LPM sup- 



pression is more pronounced at low uj, where the photon formation length (|3.25| ) is 
larger. 

In the case of usual matter, the LPM effect was observed experimentally for 
energetic electrons crossing thin metal foils, first at SLAC [42] , and more recently at 
the CERN SPS [43]. An accurate description of the data based on rigorous theoretical 
calculations is available [24,44]. 

Intermediate region X <^ L <^ L* 

The region A ^ L ^ L* is transitional between the BH and LPM regimes. The 
energy loss in this region is the same as that induced by a single effective scattering 
of typical momentum transfer qj_{N) ~ (A^lnA^) /i^, with N = L/X. We obtain from 
( P.10|) the energy loss 



AE{X <^L<^L*) ~ aElni ] , (3.33) 



which in analogy to ( p.22| ) arises from photon formation lengths 



^ =^/Ln«^/«4- (3.34) 
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The result (|3.33|) holds as long as the radiated photon "sees" only one effective 
scattering center, i.e., as long as ^/Ij^jn ^ L. This condition is precisely equivalent 
to L < L*. 

The regimes L <^ A and \ L <^ L* smoothly match when L \ (see ( p.llj ) 
and ( |3.33| )). Comparing now ( |3.33| ) and ( p.29| ), we see that the former estimate of 
AE{L = L*) involves an extra logarithmic factor and dominates at this point. It is 
associated with photons radiated outside the medium. Adding to ( ^.331 ) the LPM 



linear term (|3.29| ), we obtain 



AE(LZL*) ~ aElnf^^^ +aE^. (3.35) 

\\ m? J L* 

The linear regime sets in when the second term starts to dominate. This happens at 
the scale 

which is somewhat larger than L*. 
3.2 Energetic muon 

What happens if the ultrarelativistic particle going through the plasma is heavy, 
M ^ fi ^ eT? As already mentioned, the intensity of the BH radiation for a 
massive QED particle is suppressed by the factor ~ /i^/M^, see (^.81). For L -C A, 
we thus have 

AE{L«X) aE^^ ^ (3.37) 

The suppression factor ~ fi'^/M'^ in ( p.8|) arises from an integral over the photon 
radiation angle of the type 

e', [ de'/e' ~ e'jel , (3.38) 



which is saturated by the angles 0^ ~ = M'^/E^. 

Note now that, for a massive particle, the expression ( p.l9|) for the vacuum 
formation length is modified, 

ijiu, e) - M) ~ \ , . (3.39) 

The characteristic formation length of emitted photons is thus of order 

W - ^/^'- (3-40) 

When <^ aET, we have i^^"^^ ^ \ ^ L, i.e., the photons are mostly formed 
outside the medium. 
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Consider now the behavior of the radiative energy loss for larger lengths. There 
are two distinct cases, depending on whether is smaller or larger than avET^. 
The appearance of this scale is easy to understand from our previous study of the 
light particle case. We showed in particular that the result ( p. 26 ) - or more accurately 



( p.29| ) - in the domain L ^ L* corresponds to an electron in-medium virtuality of 
order k^^^ ~ aV ET^. We thus expect strong modifications of AE{L), as compared 
to the light particle case, when ^ avET^, and milder modifications when 
< aVET^. 

A) < aVET^ 



The condition -C ayET^ is equivalent to saying that the formation length 
( CT) exceeds the scale L* given in ( ^1231) or (^1281) , 



» L*^^— ^ « ~ aVETK (3.41) 

When L somewhat exceeds A, the medium still acts as a single effective scattering 
center, transferring the momentum q\ ~ (L/A)/i^ to the massive particle. The energy 
loss will be given by aE times the suppression factor q^/M"^ ~ (L/^) ■ (yU^/^^) as 
long as q^/M"^ -C 1. The result coincides with ( p.37| ), which is thereby valid up to 
the scale L** defined by 

L** = A— >A. (3.42) 

We can thus write: 

AE{L«Ln aE-^ L . (3.43) 

As soon as the suppression factor disappears, the problem is reduced to the 
already discussed case of light particles. The vacuum formation length ( |3.4(J| ) is 
modified in the medium to L* . For large lengths L ^ L*, the dependence of AE on 
L is linear with the same slope as for light particles, 

A^(L>L*) ~ aEL/L*. (3.44) 

As for a light particle, we can infer that in the intermediate region L** <^ L <^ L*, 
we have (see ( p.33D ): 

AE{L** « L « L*) ~ aE In ~ «E In . (3.45) 

The overall behavior is displayed by the solid line in Fig. ^. It is similar to the 
case of a light particle (Fig. with the scale L** playing the role of A. Note also 
the relation 

tp"'' , (3.46) 
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Figure 5: Radiative energy loss of a heavy charged QED particle. Solid line: moderately 
heavy particle, fi^ ~ aT^ < < aVET^. Dashed line: ~ aVET^. Dash-dotted 
line: ^ aV ET^. Dotted line: Collisional loss. 



which is analogous to (|3.23| ). 

B) > olVeT^ 

This case corresponds to the formation length ( p.40| ) being smaller than L*, 

^ M^:$>aVEf^. (3.47) 



If the vacuum formation length ( 3.401 ) is less than the would-be in-medium for- 



mation length L*, the LPM effect never plays a role. The linear law ( p.37|) is valid 



for all lengths. Indeed, it naturally extends up to the scale L = E/M"^, where the 
suppression factor (L/X) ■ (/i^/M^) ~ L/L** is still there. When L further increases, 
the suppression factor stays frozen, but the number of emitted photons grows, 

AE(L > E/M^) ^ ( aE ~ aE --^ . (3.48) 



This coincides with ( 3.37 ). The behavior of AE{L) in this case is represented by the 
dash-dotted curve in Fig. |. 

We summarize the results of this section by briefly discussing the different slopes 
which appear in Figs. § and |^. Let us consider a fixed energy E and progressively 
increase the particle mass. When the particle is light (Fig. ^ the slope for L -C A 
is larger than the slope for L L* by a factor ^jEjT . When the mass increases 
(Fig. ^, solid line), the linear regime at small L extends now to L** (instead of A) and 
has a slope reduced by a factor [i^ jM'^ . This slope will still be larger than the slope 
at L > L* as long as < a^KT^. When ~ a^ET^ (Fig. |, dashed line), 
L** coincides with L*, the two slopes also coincide and AE{L) is given by ( |3.37| ) for 
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all L. When M further increases (Fig. dash-dotted line), AE{L) is still given by 
( p.37|) , in particular its slope decreases as ~ 

Comparing this slope with the slope ~ a^T^ for collisional energy loss (see sec- 
tion 1^), we see that the radiative losses of a massive QED particle parametrically 
dominate over collisional ones (dotted line in Fig. ^ when the mass is not too large, 
or equivalently at large enough energies: 

< aET ^ E > — . (3.49) 

aT 

The latter condition for the dominance of radiative losses in a QED plasma is actually 
equivalent to the vacuum formation length i^'^'^^ E / being larger than the mean 
free path A. It is also interesting to note that (|3.49|) can be put in correspondence with 
the case of a massive particle crossing usual matter, see (|3.18|) with the replacements 
Z ^ 1 and m ^ T. 



4. Radiative loss of an "asymptotic parton" 

The experiment where an incoming energetic parton (constituent of an asymptotic 
hadron) enters a preexisting QGP, suffers radiative energy loss, and then escapes the 
QGP is, unfortunately, not feasible. We nevertheless view the problem of evaluating 
the radiative loss of such an "asymptotic parton" crossing a QGP as an instructive 
and useful theoretical exercise. We will spend some time on it before going to the 
more realistic case of a parton produced in the medium in section |^. 

The problem of the energy loss of a parton coming from infinity is somewhat 
better posed (at least, at the level of a thought experiment) for a tagged heavy 
quark than for a light parton. When the quark is heavy {i.e., when the quark mass 
satisfies M 3> Aqcd) one can think, as was discussed in the Introduction, of a 
heavy meson (or heavy baryon) scattering off a thermos bottle containing a QGP. 
The heavy quark energy loss in the QGP will be roughly the same as the energy 
difference between the initial heavy meson and the final fast hadron with the same 
flavour. 

If the projectile is a light meson (say, a pion), it contains at least two light valence 
quarks, and the energy loss through the QGP of a single quark is not observable, 
because of the absence of tagging. One can still imagine a thought experiment 
where the net energy loss of the light projectile constituents passing through some 
medium could be observed. Consider the situation where the two valence quarks of 
the incoming pion materialize as two separate jets after crossing the medium. For 
the jets to be well separated, they should have a large relative transverse momentum. 
This can happen if the pion enters the medium in a compact configuration {i.e., with 
a large relative transverse momentum between the quarks) and loses coherence due 
to in-medium rescatterings. Although the sum of the total energies of the final jets 
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source 



Figure 6: The diagram witli a tiiree-gluon vertex contributing to the third term of ([4.1j). 

should coincide with the initial pion energy, the energy loss of the light quarks after 
pion dissociation should affect the energy distribution of the leading hadrons within 
the jets. The pion diffractive dissociation process in cold nuclear matter, vr + A — 
2 jets + y4, has been studied experimentally [45], as a tool to access the pion wave 
function [46] . Light quark energy loss in nuclear matter should in principle affect this 
process. One can contemplate a similar experiment with a QGP. Thus, the energy 
loss of an "asymptotic light parton" crossing a QGP is in principle observable. 

In the following, we will successively study the cases of a light and heavy parton, 
light and heavy referring to whether the parton is lighter or heavier than the Debye 
mass yU ~ gT in the QGP. We will assume /i ^ Aqcd? so that within this definition 
a light parton, m <^ /i, can also be heavy in the usual sense, m ^ Aqcd- We follow 
the same lines as in QED (section ^ and study the behavior of the radiative loss in 
different regions of the distance L travelled by the colored particle. 

4.1 Light parton 

We will consider the case of a light quark rather than of a gluon, since it is technically 
more convenient (we will also discuss heavy quarks and it will be instructive to 
see what is changed when the mass is increased). In the following formulae, the 
parameter m stands for the quark mass when Aqcd itl <^ fi, but should be 
replaced by Aqcd if the quark is really light in the usual sense, m < Aqcd- All 
the results for a gluon are qualitatively the same as for a light quark, with the 
replacement m —>■ Aqcd and with modified color factors - but those are not our 
concern in this paper. 

i <C A: Bethe-Heitler regime 

For L ^ A, it is sufficient to determine the radiative loss induced by one scat- 
tering, as in (prrp. The gluon radiation amplitude A^rad induced by a single elastic 
scattering includes, besides the graphs of Fig. ^ (with proper color factors), also the 
graph with a three-gluon vertex shown in Fig. ^ This graph is evaluated for instance 
in Ref. [47] in the case of a massless quark. The generalization to a massive quark is 
trivial and the sum of the three diagrams gives for uj <^ E {we will work consistently 



in this approximation, as we did in QED) 



A^rad OC 



6*2 + 6*2 



e' 

0,2 02 



6" 
6'"2 + 6*2 



(4.1) 



In the latter expression, the parameters 9 = k±/uj, 6' = 6 — Og = 6 — q±/E and 
9m = m/E are the same as in ( p.5|) . We have also denoted 9" = 9 — 9g, with 
9g = q±/uj. The parameter 9g can be interpreted as the scattering angle of the 
virtual gluon of energy ~ a; in Fig. ^. The color factors can be conveniently grouped 
in terms of the commutator and anticommutator {/!:", t^} of color matrices. 

Eq. (|4.1|) can be rewritten as 



A^rad OC 



(4.2) 



where Je (already given in ( ^TSf )) and Jq read: 



J. 



Jn 
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(4.3) 
(4.4) 



The soft gluon radiation intensity can be obtained by squaring ( [4.2[ ), sum- 
ming/averaging over color indices, and normalizing by the elastic scattering cross 
section. For Nc quark colors, we obtain: 



d/rad = -,— d'9 |iv, (J, - ur + 



^/^d^9iN.J. 



1 



2N, 



72 Nc 



yJq Je) Jq 



(4.5) 



In the latter expression, the terms have been organized to facilitate the integral 
over angles. Considering a light quark with m <^ /i and using (|3.?| ), we see that the 
first term of (|4.5|) contributes to the energy spectrum as 



dL 



it) 



rad 



broad 



tts In ^ ~ In I — — 



9 9 / ' 



(4.6) 



the logarithm arising from the broad angular domain 9m 9 <^ 9g. The second term 



of ( |4.5| ) is similar to the QED case, see ( p.4|) . Its contribution to the spectrum is 



dl 



UJ 



rad 



duj 



1 1 

«s m — - ~ a, m — - 



^2 



(4.7) 



arising from the narrow angular domain 9m 9 9s- The third term of ( |4.5| ) 
does not bring any logarithm in the energy spectrum and will be neglected in the 
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following. By integrating (|4.6| ) and ( [4.7|) over cu, we see that the terms oc and 
oc contribute similarly to the radiative loss/^ 

AE-i,~a,Eln4- (4.8) 

Thus, in the light quark case m <^ /i, we have, similarly to the QED case (see (|3.1j ) 
and (CT)), 

AE(L«A) ~ yAE^:,^,, ~ a.E:^ln4- (4-9) 

To recapitulate, at fixed <^ i?, we have Og <^ 6g, and the differential energy 
spectrum receives two logarithmic contributions: a QED-like contribution, from the 
narrow region Om ^ <^ Og, and a contribution specific to QCD, from the broad 
region 9m 9 <ti 6g. The second contribution (note that it dominates at large Nc) 
reads 



a;ri£^ d'9r\ (4.10) 

broad i 

Using k± ~ co'6' and neglecting the quark mass, the spectrum ( [4.10| ) can be expressed 



as 



U!- 



duj (Pk 



^ , (4.11) 

which is the well-known Gunion-Bertsch spectrum [47]. 

We emphasize that the spectrum (|4.10|) is obtained from the QED spectrum 
( p.4|) by replacing the electron scattering angle 9^ = q^/ E ~ ^/E by the virtual 
gluon scattering angle 9g = g^/cu ~ n/oj (compare ([4.3|) and ([4.4|) ). This fact has 
interesting consequences. 

(i) At fixed u <ti E, the broad angular domain 9m ^ 9 <^ 9g contributing to ( [4.10| ) 
translates to the interval in gluon formation lengths 

4«^/(^)«^, (4.12) 

to be compared to ( ^.21 ) in the QED case. Thus, in QCD the gluon starts to 
be formed at much smaller lengths, 

(4-13) 

than the photon in QED. Since ^/('^)|jnin oc in QCD, we expect the LPM 
suppression of energy loss to increase with increasing u (in QED, the opposite 
was true, see the remark after ( p.32| )). We will discuss this in more detail below. 



-'^^This is because in the integrated spectrum we have uj ^ E, and the broad and narrow angular 
domains then coincide. 
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As long as uj ^ E in the integrated spectrum (in which case the broad QCD 
domain 9 <^ 9g coincides with the narrow QED domain 9 <C 9s)-, the typical 
formation lengths in QED and QCD do not differ (see ( p.22| )), 



4<^/(^~^)<^- (4.14) 

We thus anticipate that the average radiative loss of a light quark has a para- 
metric dependence similar to that of an electron, for all travel distances L. 

(iii) On the other hand, if the characteristic frequencies uj in the integrated spectrum 
were much smaller than E., the broad angular domain would indeed be broader 
than the narrow one, and we would obtain different parametric behaviors for 
the radiative loss in QCD and QED. We will see in the second part of this 
section that it is exactly what happens for heavy quarks. 



L > X*: LPM regime 

The derivation of the light quark energy loss when L ^ L* follows the same 
lines as in QED, with some differences which we discuss below. First, similarly to 
the photon case, for L ^ L* a single gluon is emitted in a multiple scattering process 
composed of ~ i^^'^i^)/^ individual scatterings. This amounts to exchanging N 
gluons in the t-channel. This A^-gluon state can be either color octet or color singlet 
(higher representations do not contribute to the quark scattering amplitude). In 
the color singlet case, the physics is the same as in QED and the radiated gluon 
is emitted within a narrow cone^^ of angle 9'^{N) ~ Nfi'^/E'^. For the color octet, 
the physics is the same as for the process with one gluon exchange discussed above 
where two cones are present, but with the narrow and broad cones angles being 
now of order 9g{N) ~ N^^ /E"^ and ^^(iV) ~ A^/i^/cu^, respectively. In the large Nc 
limit, the probability to have a singlet t-channel A^-gluon state is suppressed by l/A"^ 
compared to the probability to have an octet. This suppression is of the same order 
as the suppression of the Abelian contribution yielding the narrow radiation cone 
in single gluon exchange (the second term in ( ^.5| )). In other words, the dynamics 
of gluon emission in a multiple scattering process is roughly the same as for single 
scattering, the only difference being that the characteristic momentum transfer /i^ is 
multiplied by the factor A^ ~ ^Y'^{uJ)/\. 

The in-medium gluon formation length P^'^^i'^) is obtained from ( p.l9|) by re- 
placing the gluon radiation angle 9"^ by 9g{N) ~ Nfi'^/u'^, 
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In contrast to what happens in the BH regime, this cone does not produce any logarithmic 



factor, see the footnote related to (3.26) 
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The gluon formation length increases with a;, an opposite behavior compared to 
QED, see ( p.25|) . The gluons with uj ^ E are formed within the length L* given in 



( |3.23| ), notwithstanding, leading to the same parametric dependence of the radiative 
loss as in QED, 



AE(L>L*) ~ ~ "«^VX^ ~ allVET^. (4.16) 

There are, however, differences between the QED case of an electron and the 
QCD case of a light quark. 

(i) Similarly to QED, due to the long tail in the Coulomb scattering potential, 
the typical transverse momentum exchange q^ypiN) after scatterings is not 
exactly N/i'^. As shown in Appendix A, it is ~ (A^ In A^) jj^ in QED, but in QCD 
this simple dependence is replaced by ( |A.25|) due to the running of as- For very 
large A^, QtypiN) is given by the expression (|A.26|) that does not involve the 
factor In A^. This means that the factor ~ \/\nE in the energy loss of energetic 
light partons in the large L region disappears for asymptotically large energies 
ln{E/T) > ln(/i/AQCD)- We have instead of (^^) and ( F^ , 



\E 



/i2 1n-^ ' 

Aqcd 



(4.17) 



AE(L>L*) ~ asE- aiL JET^ \n-^ . (4.18) 

L* Y Aqcd 

When the energy is large but not asymptotically large, the dependence is more 
complicated. 



AE{L^L*) ^ alL 



(ii) Due to the difference between the gluon ( [4.15| ) and photon ( p.25| ) formation 



lengths, the LPM gluon energy spectrum reads, instead of ( pj.30D 



where ujc is still given by (p.31|) , but depends now on the parameters /i and A 
appropriate to a non-Abelian plasma. As expected, the LPM suppression in 
QCD increases with increasing cu, contrary to QED. Note that below AyU^, the 
formation length is less than A and we are in the incoherent BH regime where 
the spectrum is given by (O) times the number of rescatterings L/\. 
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Inserting the logarithmic factor discussed above, we obtain 



^^(L) ~ alLJ^ln-J^ , (4.21) 
duj y Aqcd 

where we assumed ln(co'/T) ^ ln(/i/AQCD) (otherwise the expression is more 
comphcated, similarly to ( |4.19| )). 



Intermediate region X <^ L <^ L* 

In this region, the reasoning used in the case of an electron can be directly 
transposed to the QCD case of a light quark, yielding the result 



A^(A < L < L*) ~ asE In ( ) . (4.22) 



As in the QED case, the result ( |4.22| ) corresponds to the medium acting as a 



single effective scattering center. The medium size dependence appears only through 
the total momentum transfer ~ (A^lniV) /x^, where N = L/X. 

Eq. ( [4.22| ) represents the so-called factorization term mentioned previously in 
Refs. [12,13,15,41]. Although the main goal of Ref. [15] is to address the radiative 
loss of a quark produced in the plasma, it is mentioned there that in the case of 
an "asymptotic quark" entering the medium, the factorization term can be dropped 
when calculating the induced energy loss. As a consequence, the result of Ref. [15] 
for AE in the region A ^ L ^ L* is AE oc agUc, instead of ( [4.22| ). But in fact, for 
an asymptotic quark there is no distinction between the induced and total radiative 
energy loss,^° and the factorization term ( [4. 22] ) should be kept. Despite the fact 
that the factorization term has a smooth (logarithmic) dependence on L, it actually 
dominates over the term calculated in [15]. Indeed, for L <C L*, we have a^tUc -C asE. 

To illustrate this point, we represent in Fig. ^ the gluon radiation spectrum in 
the region \ <^ L <^ L*. For L <^ L* we have Uc ~ L^yU^/A <^ E, and the spectrum 
is given by (|4.20 ) or ( 4.21 ) in the interval A/i^ <^ ^ cjc- As mentioned above, for 



UJ <^ Xfi^, the spectrum is given by ( [4.6| ) times L/X. For uj ^ Uc, the formation time 
of the gluon is larger than L, and the spectrum is the same as for a single effective 
scattering, i.e., it is obtained, again, from ( |4.(j| ) by replacing /i^ fi'^L/X. This flat 
part of the spectrum gives the dominant contribution to the light quark energy loss. 



see ( ^1221) . 

Thus, the law AE{X ^ L ^ L*) oc [13-15, 17] for the induced energy loss of 
a light quark produced in the medium (we will review this case in Section is not 



^"^This is because an "asymptotic, on-shell quark" does not radiate in the absence of the medium. 
This is in contrast with a quark produced in a hard subprocess, which radiates even in vacuum. 



See section |^ for more discussion on this point. 
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Figure 7: Gluon radiation spectrum by an asymptotic light quark crossing a hot QCD 
medium of size L, with X <^ L <^ L* (double logarithmic plot). 



valid for an asymptotic quark. Combining ( ^l9|) , ( [4.22| ) and ( [4.16|) , we see that the 



average radiative loss of an asymptotic light quark is similar to that of an asymptotic 
electron crossing a QED plasma. It is represented in Fig. ^. 

4.2 Heavy quark 

We have found that the radiative energy loss of an asymptotic light quark crossing 
a QGP has the same parametric dependence (apart from a logarithmic factor when 
L ^ L*) as for an electron crossing a hot e~^e~ plasma, despite drastically different 
radiation spectra in these two cases. We will see now that the heavy quark radiative 
loss is different from the heavy muon loss in QED, and has a richer parametric 
dependence. The non-Abelian dynamics manifests itself more clearly for heavy than 
for light quark radiative energy loss. 

Consider first the BH regime, L <^ A. 

For heavy quarks, M yU, the gluon radiation intensity is suppressed, similarly 
to heavy leptons in QED. However, in the non-Abelian case, the suppression (dead 
cone effect) is not so strong since soft gluons are emitted in a cone broader than in 
QED, see ( [4.1CI| ) and the related discussion. The estimate ( |3.9| ) for the QED radiation 
spectrum should be replaced by 



dliad 

UJ 



duj 



^~a.ln(l + ^^- -^^^ ■ ^'"^ 

Iscat. -J V 



as I rf^^J/~a.ln(l + ^l ~a.ln(l + ^— ). (4.23) 



Integrating the latter spectrum over u and multiplying by the scattering probability 
L/Xwe obtain 
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Thus, for L <^ A the radiative loss of an asymptotic heavy quark is suppressed by a 
factor ~ /i/M, instead of ~ /i^/M^ in the QED case, see ( p.371 ). A very important 
distinction of QCD compared to QED is that the spectrum ( |4.23| ) of a heavy quark 
is soft. Indeed, the characteristic energy of emitted gluons is eschar ~ fJ,E/M -C E, 
and the spectrum falls rapidly (as l/uj"^) beyond this scale. (Instead, the photon 
radiation spectrum of a heavy charged particle remains flat: mass effects bring about 
a uniform suppression ~ /x^/M^ for all u.) The result ( ^.24| ) arises from small gluon 
energies, u ~ cUchar- The typical gluon angles contributing to (|4.24| ) are of order 



^char ~ ~ /^^/^char ~ ^M" Heucc the characteristic gluon formation length is 

tfQCD heavy 1 ^chav E 



(4.25) 



Similarly to the Abelian case, the behavior of the radiative loss at larger lengths 
depends on the ordering of different length scales: the mean free path A, the charac- 
teristic formation length, or the scale L** where the suppression ~ /x/M disappears. 
The ordering of these scales depends on the heavy quark mass. In the Abelian case, 
we had two mass regions ( |3.41| ), ( ^.47] ). In the non- Abelian case, there are three 
distinct regions. 

A) < ctsVET^ 



In this case, the smallest length scale is A ~ l/{asT). For L -C A, the law ( [4.24| ) 



holds. When L exceeds A, the energy loss will be the same as that induced by one 
effective scattering of momentum transfer /i^g ~ (L/A)/x^, 

L ~ asE . 4.26 

^ ^ M M ^ ' 

This law is valid in the region 

A « L « L** = A— ~ ^ . (4.27) 

At the scale L**, the suppression ~ [le&jM disappears and the physics becomes the 
same as for light quarks. When L exceeds L**, but is still less than the in-medium 
formation length L* given by ([4.17|) , there is still one effective scattering, and we 



are in the intermediate region. When L ^ L*, the number of effective scatterings 



grows as L/L*, and the radiative loss is given by the light quark estimate ( [4.18|) . A 
schematic plot of the energy dependence is drawn in Fig. R. 



B) (XsVeT^ < < ctsE^ 



When exceeds OsV ET^, the scale L* becomes smaller than the scale L** 
where the dead cone suppression disappears. Before proceeding further, note that the 
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Figure 8: Radiative energy loss of an asymptotic heavy quark with <C Og V ET^ . 
{Dotted line: colHsional loss.) 

estimate L* ~ ^JXEJ^ for the in-medium formation length does not hold anymore 
in this case, and L* should be replaced by another scale L* . 

Indeed, recall that the in-medium formation length is defined by the condition 
that it coincides with the vacuum formation length for an effective scattering of 
transfer /i^g ~ (L/A) /i^. For a light quark this condition reads 

L ~ ^ ^ L* ~ J^. (4.28) 

/^eff V ^ 



For a massive quark the vacuum formation length is given by ( |4.25| ) and we obtain 



E ~ ( \E^V'^ fU^V''' If E V'^ , , 



Note that the condition agVET^ ^ ^ defining the mass interval under 

consideration is equivalent to A -C L* <^ L*. 

For L ^ L*, the size of an effective scattering center responsible for the emission 
of one gluon is ~ L*, and the dead cone suppression factor ~ fies/M stays frozen at 
the value 

r~ / — \ 

/AteflF\ ^ JJ_ ^ I asVET^ 



\M J max MV A \ M2 



(4.30) 



The energy loss displays in this region a linear dependence on L with the slope [48] 

2/3 / X 2/3 



AE{L » t) ~ a,E i- (^) ~ a,L ( ^) ~ T'L ( ^) 

^ ' L* VM/max \XM J ' \MJ 



(4.31) 
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Figure 9: Gluon radiation spectrum of a heavy quark produced in a hot QGP for 

asVET^ « < asE"^ and for L-^L*. 

The schematic plot of AE{L) looks as in Fig. ^, but with two qualitative dis- 
tinctions: (i) the scale L* is replaced by L*; (ii) in the region between A and L*, 
AE{L) oc VI (as in dOep ) instead of AE{L) oc hiL as in (|33^) . 



Let us discuss now the energy spectrum of emitted gluons. Consider L ^ L* . 
We expect the spectrum to differ from the light quark case (see ( [4.20| ) and Fig. 0) 



as soon as the typical angle Otyp contributing to ( [4.2CI| ) becomes smaller than the 
parameter 6m = M/E. The typical angles associated with (^4.2CI|) read 



where we used ( [4.15| ). The condition 6'typ < 9m is thus equivalent to 



Li^/A /,,2z?4\l/3 z?2\2/3 

< — >Q 
cu-^ E'^ 

Note that L ^ L* is equivalent to a) ^ Uc- When u > to the spectrum is given by 



dl L n^E^ 1 , 

uj—— ~ a^^^ — A ■ ~ Qs-L , , • ^ ~ Cis n — [UlI>UJ) . (4.34) 

duj ^ uj^ ^ ' ^ ' 

For UJ < UJ the spectrum is given by ([4.20|) . The overall spectrum is depicted in 
Fig. p. Note that the behavior oc uj~'^ that we find for uj > uj differs from the 
behavior oc uj~'^l'^ obtained in Ref. [26]. Integrating the spectrum we recover the 
average loss ([4.31|) , the latter being dominated by ~ cD. Noting that 




Cj ^ E ~ — \/ — , (4.35) 
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the average loss ( 4.31 ) can be rewritten in a compact form, 



AE(L>L*) -a.cl^i, (4.36) 

L* 



to be compared with the hght quark estimate (|4.16|) . 



C) > (x.E"^ 

Finally, for very large masses, the LPM effect does not play any role and the BH 
linear law ( |4.24| ) holds for a// lengths. Indeed, when exceeds UgE'^, the formation 



length L* becomes smaller than A and gluons are emitted in individual incoherent 
scatterings. In this case, radiative losses are suppressed compared to coUisional ones, 
which can be checked by comparing the slope in (f4.24|) with the slope ~ alT"^ for 
coUisional loss. Because of this, we stress here that the radiative loss ( |4.24| ) might 
be difficult or impossible to observe (even in a thought experiment). See section ^ 
for more discussion of this point. 

As was mentioned above, the spectrum of emitted gluons is soft in this case and 
given by ( [4.23| ) multiplied by the number of scatterings L/\. The spectrum starts to 
fall down as ~ cu"^ at the scale nE/M rather than at cD, as it did in the intermediate 
mass region. 

The linear energy dependence AE{L) oc E may suggest the description of heavy 
quark radiative losses in terms of radiation length, as it is usually done for ultrarela- 
tivistic electrons in usual matter. But it is not convenient here for two reasons: (i) In 
contrast to ultrarelativistic electrons, here coUisional losses dominate over radiative 
ones; (ii) The spectrum of emitted gluons is soft and energy loss fluctuations are 
much smaller than e.g. for electrons forming atmospheric showers. 



5. Radiative loss of a particle produced in a plasma 

Here we consider the case of a charged particle produced inside a plasma. This 
situation is much more natural for QCD, where an energetic parton can be produced 
in a hard partonic subprocess inside the hot medium formed in heavy ion collisions. 
We however start by considering the less natural but simpler QED case of an electric 
charge produced in a QED plasma. We will then study the QCD case. 

5.1 Hot QED plasma 
5.1.1 Electron 

When an energetic charged particle is created in a hard process, it radiates brems- 
strahlung photons. This radiation occurs even when the particle is created in vacuum, 
and should be distinguished from the medium-induced radiative loss. 

In this section, we will consider the situation of a fast and light charged particle 
produced in a hot QED plasma (think for instance of deep inelastic scattering off a 
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QED plasma, or of direct lepton production in a QGP), and focus on its medium- 
induced radiative loss. This loss is associated with the components of the particle's 
radiation field coat which had the chance to be formed within the distance L travelled 
by the newborn particle in the plasma, such that those components can be released as 
emitted photons during subsequent rescatterings. In other words, only the photons 
whose formation length ( ^.KJ] ) does not exceed L contribute to the induced radiative 
loss, 

1 

In the case of an asymptotic particle studied in sections ^ and ^, we found some 
contributions to the energy loss arising from if ^ L. In particular, in the BH 
regime L -C A, the electron energy loss (|3.11| ) arises from photon formation lengths 
if ^ E / ij? ^ L, see ( |3.22| ). Due to the prescription ( [5.1| ), this contribution should 



ff{uj,e)^ — <L. (5.1) 



now be disregarded. We should only count the photons whose formation length does 
not exceed L, which brings an additional suppression in medium-induced radiative 
losses. Thus, when a light particle is created inside the plasma, there is no BH regime 
whatsoever. We will shortly see into what kind of behavior it is transformed. 

On the other hand, the result ( |3.26| ) (or, more accurately, ( |3.29| )) for the radiative 
energy loss in a thick medium, L ^ L*, should also be valid when the particle is 
created in the plasma rather than in the remote past. Indeed, the result (|3.26|) arises 



from photon formation lengths i^'^'^i^) ~ i^'^'^iE) ^ L* <^ L, thus satisfying the 
prescription (|5.1| ). When L ^ L*, the particle forgets the conditions of its birth. 

Consider first the region L <^ \. In this case, the particle undergoes one scat- 
tering with probability ~ L/X. The photon emission amplitude is given by the sum 
of the two diagrams in Fig. For small photon frequencies uj <^ E and small 
scattering and emission angles, it can be evaluated as oc ee ■ J{L) with 

= • (5-2) 

Here L is the distance travelled by the particle between its production and scattering. 
We assumed that the particle is massless. The term oc 6' is the contribution of the 



graph in Fig. 10a and the term oc 9' is the contribution of the graph in Fig. 10 b. 



The result ( p.2| ) for the amplitude is rigorously derived in Appendix B, but its 



structure looks rather natural in the context of the above heuristic reasoning. When 
L is large compared to the formation length (|5.1| ), one can drop the rapidly oscillating 
factor ~ exp{—iujL9'^/2} and the current J{L) is reduced to the expression ( p.5|) for 
an asymptotic particle. On the other hand, for a photon formation length larger 



than L, the contribution of the graph in Fig. 10a is suppressed and only the second 



graph remains, corresponding to photon emission from the final electron line - as in 
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For heavy particles, this is not so, see the discussion below. 
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Figure 10: Photon radiation of an electron produced and scattered in a QED plasma. 



the absence of rescattering. Since we are interested in the medium-induced radiation 
intensity, we should subtract the latter contribution (at the cross section level), giving 



UJ- 



dl 
I — 

duj 



induced 
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The averaring is done over the transverse momentum gl exchanged in the scattering 
(we remind that 9' = 6 — 6s = 6 — q^/E). Let us average first over the azimuthal 
directions of 9s. Using the identity 
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we obtain 
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We should now average over 9s with the normalized probability P{9s), obtained 
by squaring the (momentum space) screened Coulomb (Yukawa) potential, see ( p.2| ). 



The result is 
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AE2 Jq 9^{9^ + iiyE^] 



(1 - cos(cuL0V2)) • (5.7) 



induced 

When ujLfi'^/E'^ <^ 1, which is true in the region L -C A we are now considering, the 
integral is saturated by the emission angles^^ 



^^The last inequality in ( |5.8| ) follows from E ^ \fp ^ T. We also assume EL '3> 1, i.e., even 
though L is smaller than A, it is still larger than the wavelength of the energetic particle. 




Figure 11: Photon energies and emission angles satisfying if{uj,9) < L. 



The integral in (|5.71) is of order ~ uL. Thus, the spectrum is evaluated as 



dl 

duj 
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induced 



A 



UJ UJUJc 

~ a- 



and the energy loss reads 



(5.9) 



(5.10) 



It is interesting to mention the analogy with the discussion of the energy loss of 
an asymptotic and massive particle in Section p.2| . The effective cutoff 9"^ > 1/ {EL) 



which arises here is similar to the "dead cone" cutoff 9"^ > M'^/E'^ in the integral 



( p.38| ), with the parameter EjL playing the role of mass. One can actually obtain 
the estimate ( ^lOD by replacing M"^ ^ E/L in ( ^^431) . 

We see that /S.E{L) has a quadratic rather than linear dependence at small L. 
The reason for that is quite transparent, as we qualitatively explained before. When 
a hard electron is just born, it has not grown its radiation field coat yet and is not 
able to radiate. Roughly speaking, the capacity dE/dx to radiate vanishes at L = 
and then grows linearly with L. The integration 



AE 







^ , dE 
dx — — 

dx 



(5.11) 



yields another factor ~ L. 

It is instructive to discuss a more heuristic derivation of the estimate ( ^.10| ) 
which does not use the exact expression ( ^.2|) of the radiation ampUtude, but simply 
consists in integrating the spectrum ( p.6|) (derived for an asymptotic particle) over 
u and 9 with the constraint 1/{uj9'^) < L. The corresponding integration domain is 
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depicted in Fig. 11. Since 6 ^ ^ /E\ the angular spectrum can be approximated 



by / The energy loss then reads 
AE(L<A) ~ a- / du de^-^ ~ / duuL ~ au, , (5.12) 

^ Jo Jl/(wL) Jo 

arising from the typical values 9^ ~ 6'^;^ ~ 1/{EL) and uj ^ E. This coincides with 



We want to emphasize, however, that although the latter argument is very simple, 
physically transparent, and gives the correct result, it is heuristic and does not repro- 
duce a certain subtle dynamical feature which is displayed in the accurate derivation 
based on (|5.2| ).^^ Indeed, the argument leading to the heuristic estimate ( |5.12|) implic- 
itly assumes that the momentum transfer is fixed at the value q\ = fi^, i.e., it refers 
to an hypothetical model where the Coulomb probability density ( |5.6| ) is replaced by 
P(6'^) = 6(9'^ — jji^ / E"^) (corresponding to a scattering potential V{r) ~ Jo{fj,r)). On 
the other hand, substituting 9"^ — > fi^/E"^, the expression ( p.5|) which follows from a 
diagrammatic analysis would yield 



dl 

duj 



which leads to an energy loss suppressed compared to the expression ( |5.12|) . Thus, 
the estimate ( ^.12|) fails for the hypothetical model where P{9l) = 5{9l — fi'^/E'^). 



For the more realistic Yukawa potential, the results based on (|5.5| ) and ( ^.12|) are 
the same, but the integrals are saturated in different kinematical regions. With the 
heuristic prescription ( |5.12|) , the characteristic radiation angle 9^^^ ~ 1/{ujL) is much 



larger than the characteristic scattering angle 9"^^.^^ ~ fi'^/E'^. In the more accurate 
formula ([5.7|), they are both large (see (p.4[)), 

~ Cat ~ ^ • (5.14) 

A distinctive feature of the Yukawa scattering potential is that it allows very large 
transfers compared to the typical transfer ~ /i. Thus, the constraints 9'^ > 9'^ (see 
( ^.4|) ) and 6*^ ~ 1/{ujL) ^ IJ^ / E"^ (see ( ^.8|) ) can be realized simultaneously. 

Let us now consider a larger medium, L ^ X, and focus on the Yukawa scattering 
potential. The electron is now scattered times, with N ^ L/X. If L ^ L*, only 
one photon is emitted. The amplitude of photon emission in the multiple scattering 
process is evaluated accurately in Appendix B. It happens that the main contribution 



^^In other words, the diagrams in Fig. |l^ describing photon emission during a single scattering 
of a particle produced in a plasma dictate a slightly different dynamics, compared to the single 



scattering diagrams for an asymptotic particle in Fig. ^supplemented by the constraint {^A). On 
the other hand, the physical arguments that derive the multiple scattering dynamics on the basis 
of the single scattering diagrams and the notion of formation length seem to work in all cases. 
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Figure 12: Radiative energy loss of an electron produced in a QED plasma. The quadratic 
dependence ( 5.15 ) at L ^ L* is replaced by the linear dependence ( 3.29| ) at L ^ L*. 
{Dotted line: collisional loss.) 



to the radiation spectrum arises from the region where one of the scattering angles 
is large as in ( p.l4| ), while all other scattering angles are relatively small ~ fi/E. 



In other words, one of the N scattering momenta is ~ EjL (assuming u E) 
and is much larger than the characteristic momentum transfer /Xefr ~ /i\/iV due to 
all other scatterings.^^ Any of the scatterings can be distinguished in this way, 
which gives the factor ~ L/A in the radiation spectrum. The latter is still given 
by the integral ( ^.7|) , with the factor L/X interpreted as the characteristic number of 
scatterings instead of the scattering probability. The estimate for the energy loss is 
still given by ( p.lO| ), which thus holds in the extended range L <^ L*, 

AE(L<t:L*) ~ auj^^ a^T^L^ . (5.15) 



Finally, for L ^ L*, the restriction (|5.1| ) is not effective, as we already mentioned. 
We obtain the same dependence ( p.29| ) as for a particle coming from infinity. The 
results are represented in Fig. |T2|. For large lengths, AE is linear in L, with dE/dx oc 
\fE which is familiar from the previous sections. For L smaller than L*, AE{L) 
displays the quadratic dependence (]5.15|). 



5.1.2 Muon 

Consider now the radiative losses of a massive particle created in the plasma. As was 
also the case for a particle coming from infinity, the behavior of AE{L) is different 
in the regions ( |3.41| ) and ( |3.47| ). 



A) < aVET^ 



24 



The condition y/E/L ^ fiy/L/X is equivalent to L ^ L*. 
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We will show that there is no difference with the case of light particles in this 



region, and the behavior of AE{L) is the same as in Fig. |T^. Recall that for a massive 
particle the vacuum formation length is given by ( ^.391 ). 



When L <^ L*, we showed previously that the induced radiative loss of a light 
QED particle arises from u E and 6*^ ~ 1/{ujL) ~ 1/ (EL). Thus, as long as 

1/{EL) » 91 = M'/E' , (5.16) 

the result ( p.lOl ) for a light particle will also apply to a heavy particle. Intuitively, 
this happens when M is small compared to the effective "mass" ~ ./eJL of the 
light particle (see our discussion below ( |5.1CI| )). To see that the condition ( |5.16| ) is 



satisfied, note that < aVET^ implies L* < E/M^. When L < L*, we have 
<^ E/L* <^ E/L. Thus, the result ( |5.15D is valid also for a moderately massive 
particle. 

The effects due to a nonzero mass are irrelevant in the range 
also for large lengths, L ^ L*. Indeed, when L ^ L*, the electron radiative loss 
arises from uj ^ E and photon formation lengths 1/{ijj9'^) of order L*, implying 
~ I /{EL*) > elj, where we used again L* < E/M"^. Thus, the light particle 
result (|09| ) is valid in the region < aVET^. 

An equivalent way to understand this is as follows. The length above which the 
mass can be neglected is determined by /igg ~ (L/X) /i^ ~ M^, i.e., by the scale L**, 
see ( [4.27| ). For <^ aV ET^, we have L** <^ L*, and mass effects can a fortiori be 
neglected when L ^ L*. 

B) > olVET^ 

In this region, the characteristic vacuum formation length ^/ ~ E /M'^ is smaller 
than the scale L* and shows up first. The quadratic law (|5.15|) extends only up to 
the scale L ~ E/A/P, after which it is replaced by the law ( |3.37| ), the same as for 
asymptotic and heavy QED particles. 

5.2 Quark gluon plasma 
5.2.1 Light parton 

An energetic, high p± light parton produced in a proton-proton collision can be 
"observed" via the jet of hadrons that it produces. Those hadrons are the products 
of the parton's bremsstrahlung induced by its sudden acceleration at the moment 
of its creation. In the process of building its asymptotic (t — >• +oo) field coat, the 
initially "bare" parton radiates quasi-collinear DGLAP gluons. The parton energy 
at the time of its production can in principle be determined by measuring the total 
jet energy. 

If the parton is produced in a finite size plasma, the parton energy loss due to 
its rescatterings in the hot medium can affect the structure of the hadron jet. In 
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particular, such medium-induced energy loss leads to the suppression of large p±_ 
hadrons (jet-quenching) in ultrarelativisitic heavy- ion collisions [9], when compared 
to proton-proton collisions. Also, medium-induced gluon radiation will enhance the 
hadron multiplicity within the jet. Here we want to derive the medium-induced 



radiative energy loss of a light quark. As explained in section |5.1| , the formation 
length of the medium-induced gluon radiation must be smaller than L, see ( |5.1| ). 

Consider first the region L <^ X. The physics is the same as in the Abelian case 
with the only difference that the characteristic Abelian radiation cone width ^ fi/E 
should be replaced by the no n- Abelian one ^ fi/uj (see section ^TT] and our comments 
below (|4.11|) ). Thus, the QED expression (|5.71 ) is transformed into 

QCD 



CO 



duj 



induced 



dO' 



(1 - cos(cuL0V2)) • (5.17) 



The spectrum ( |5.17| ) has two different forms depending on whether the gluon vacuum 
formation length (!.f{uj) ^ uj/ ^'^ is smaller or larger than L, 
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The contributions of these regions to the induced energy loss are 



A£;qcd,i(^ < A) 



A-Eqcd,2(-^ < A) ~ a^tUcln 



E 

Li? 



(5.18) 
(5.19) 

(5.20) 
(5.21) 



The second contribution is logarithmically enhanced. 

The expression ( |5.21| ) for the light parton radiative loss in the BH region was 
obtained in Ref. [19]. Its origin is the same as in QED. Namely, the energy loss arises 
from radiation angles larger than the typical scattering angle, 6"^ ~ 1/{ojL) ^ f^'^/E'^ 
in QED, and 0^ ~ 1/ (ooL) ^ /^^/t^^ in QCD. Thus, as was also the case for the QED 
expression (|5.1CI|) , the QCD loss (|5.21|) is specific to a Yukawa scattering potential. 

In the region \ <^ L <^ L* and for small enough frequencies, the medium effects 
come into play and the formation length is given by ( |4.15|) . One can distinguish two 
(or two and a half, if you will) regions in the spectrum. 



i'^^^{oj) <^ L, implying oj <^ oJc, the spectrum is the same as for the asymp- 



If 

totic particle in the LPM regime, see ( 4.20 ) and Fig. |^, 
dl 



OJ 
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dl 
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induced 
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induced 
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[OJ < A/i ) 
(A/i^ < OJ < OJc) 



(5.22) 
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The latter arises from typical emission angles as given in ( [4.32 ), 



1 .MA 1 



^ ^;;;;:dr^ ~ V^-V > (5-23) 



uif'^'^iuj) V ujL 



The contribution of the region u < X^^ to the energy loss is small. The region 
Xfi"^ < u < Uc yields 

A^i(A<L<L*) ~ asuj^. (5.24) 

This contribution is specific to QCD (in QED, i^'^'^i^) exceeds L if the latter 
is smaller than L*) and was identified in Refs. [13, 17] 



25 



(ii) When u > uJc (note that tUc -C -E as long as L <^ L*), if'^lu) > L. In this case, 
the radiation spectrum is the same as for a single effective scattering of typical 
scattering angle N^x^ /uP, with N L/X. The spectrum is then estimated from 
(|5.17|) . The estimate (|5.19|) follows, but is now valid when l/^uL) ^ N^'^/uj^ 
(which exactly coincides with the condition u ^ Uc), 



dl 

doj 



~ a,— {uo^<uj <E) . (5.25) 



induced ^ 



This part of the spectrum arises from the typical angles 

^ (5.26) 

ojL 

(meaning that the formation length is of order L) , and contributes to the energy 
loss as 

E L* 
Ai?2(A < 1/ < L*) ~ a^cjcln — ~ a^cjcln — . (5.27) 

LUc L 

This contribution was discussed in Ref. [19]. It is the QCD analog of the QED 
expression (|5.15|) and depends on the presence of a long high-momentum tail 



in the Coulomb scattering potential. 



The full radiation spectrum is represented in Fig. |13|. In the region uj <^ cjc, it is 
the same as for an asymptotic quark (see Fig. In the region uj ^ Uc, the spectrum 
is suppressed compared to the case of an asymptotic particle due to the constraint 
> l/(ci;L). This region still gives the dominant contribution to the energy loss.^^ 
At large energies, the contribution ( p.27| ) dominates over that of ( |5.24| ).^'^ 



^^In (5.24) we have not displayed the logarithmic factor ~ ln{L/X) which might be associated 
with this contribution. This factor arises in Refs. [13, 17] for a Coulomb scattering potential. But 
we cannot say with certainty that its presence is a model-independent statement. 

Note that other quantities than the average energy loss, involving all moments of the distribution 



26 



dl/duj, such as quenching factors [49], might be dominated by the soft part ( 5.22 ) of the spectrum. 
^^This holds irrespectively of whether (5.24) involves a factor ^ ln(iy/A) or not. 
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Figure 13: Induced gluon radiation spectrum of a light quark produced in a hot QCD 
medium for X <^ L <^ L* (double logarithmic plot). 

Finally, when L ^ L*, the parameter ujc becomes larger than E and ceases to 
play a role. The spectrum is given by (|5.22D for all u > A/x^, as for an asymptotic 



particle. Integrating over u reproduces the estimate ( [4. 161) for the asymptotic particle 
energy loss in a thick medium. 

The overall behavior of AE{L) is the same as in the Abelian case,^^ see Fig. 
Let us stress again that the quadratic growth at small L ^ L* is a specific feature 
of the medium-induced energy loss of a newborn particle. It displays itself both in 
Abelian and non-Abelian plasmas. 

5.2.2 Heavy quark 

Finally, let us discuss the radiative losses of a heavy quark created in the plasma. As 
was the case for an asymptotic quark, there are three main mass regions. However, 
the case of a heavy quark produced in the plasma is more complicated. As we 
will see, the second mass region, namely agV ET^ -C <^ asE'^, splits into two 
sub-domains, where the logarithmic dependence of AE is slightly different. 

A) < asVET^ 

In this region, the dependence of AE{L) is the same as for light quarks. The 
reason for that is basically the same as in the Abelian case. The characteristic gluon 
radiation angle, which is of the same order as the photon radiation angle, as in ( |5.26| ), 
or exceeds it, as in (|5.23| ), is much larger than 9\j in the whole range of L and 

B) (XsVeT^ < < 



UJ. 



28 



As far as AE{L) is concerned, the only difference between QED and QCD is the smooth 



logarithmic enhancement contained in (5.21) and (5.27) 
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In this range, the relevant in-medium formation length L* given in ( [4.29 ) is larger 



than A, but smaller than the light quark in- medium formation length L*. The law 
AE{L) ~ valid at small L is replaced by the linear dependence at the scale L* 
rather than L* . This is the main effect brought about by the quark mass. 

In addition, there is a more subtle effect: the logarithmic factor multiplying 
agUOc in the estimate of AE{L) at small L might change. Indeed, the logarithmic 
factors in (5.21) and ( |5.27| ) come from integrating the spectra ( p. 191 ) and ( [5. 251 ) over 



the intervals L^i^ < uj < E and Uc < uj < E, respectively. But a (large enough) 
mass brings about an effective cutoff in the spectrum when the characteristic angle 
~ 1/{ujL) becomes smaller than 6lj. This happens when the gluon energy exceeds 
the scale 

- - (5.28) 

Beyond the scale a;^, the spectrum falls down rapidly as ~ l/u^. This is a character- 
istic behavior of the hard part of energy spectra beyond the mass-induced cutoff, cf. 
( [4.23|) and Fig. ^ Evidently, the statement above makes sense only when uJn < E, 
i.e., when L > E/M'^. The results for AE{L) are thus slightly different depend- 
ing on whether E/M^ > A {i.e., < a^ET) or E/M^ < A {i.e., > a^ET). 
Representing the radiative loss as 

AE{L <^ L*) asUclnR , (5.29) 

we quote the estimates for R in the two relevant subregions.^^ 



Bl) (XsVET^ < < otsET 

In this subregion we have 



R 



Bl 



r E/{L^^^) {L « A) 

I E/uJc (A<L<E/M2) . (5.30) 
[ LUa/uJc {E/M^ <^L<^L*) 



For illustration, we represent in Fig. |T1| the induced radiation spectrum in the last 
case, namely E/M"^ < L < L*. 

B2) a^ET < < a^E^ 

Here the estimates for R become 

f E/Lfi^ {L<t:E/AP^ 



R 



B2 



\ cun/( V) (^/^2 ^L<^\) . (5.31) 
\ unloJc (A < L < L*) 



^^We do it for completeness, although the modification of the logarithm's argument is probably 
a too subtle effect to be observed in experiment. In addition, there might be logarithmic factors 



not coming from J dLo/uo, see the footnote after Eq. ( ^.24 ) 
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Figure 14: Induced gluon radiation spectrum of a heavy quark produced in a hot QGP 
for UsVET^ «; M2 < asET and E/M^ L <^ L* (double logarithmic plot). 

The logarithmic enhancement in AE{L) disappears at L ~ L*, where the energy 
loss is 



2/3 



AE{L ~ L*) ~ as uJc{L*) ~ a^E ( ) , (5.32) 



the same estimate as ( [4.30| ), ( ^.31| ) for an asymptotic heavy quark. As was mentioned, 



when L S> I/*, the quadratic law ( |5.29| ) is replaced by the linear one with the slope 
given in ( |4.31| ). 



C) > ctsE^ 

When the mass is so large, the scale L* becomes smaller than A. In this case, 
medium effects do not affect the formation length and the latter is given by E/{^M) 
(see ( [4.25| )) rather than E/^fiesM) (see ( [4.29D ). At this scale, the quadratic law in 



L is replaced by the linear law ( [4. 241) , the same as for an asymptotic heavy quark. 

At L -C E/{fj,M), the energy loss is estimated as in ( p.29| ). As in the previous 
case, the argument R of the logarithm depends on whether L < E/M"^ (in this case, 
the light quark estimate (|5.21|) still holds) or L > E/M"^ (in this case, the upper 



cutoff ( I5.28D is introduced in the spectrum). To recapitulate. 



2^ 



r E/( V) (L « E/M , 

^ \uJul{L^^) {E/M^ <^L^E/{^iM)) • ^ ' 

We observe that the law AE{L) ~ asCUc oc is universal and is not modified 
at small enough L, however large the quark mass is. On the other hand, the larger 
the mass, the earlier the change of regime between the quadratic law and a linear 
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Figure 15: Induced radiative energy loss of a heavy quark produced in a QGP. A) <C 
B 



OsVET^ <C <c q^£;2. ^2 ^ Q,^^2^ {Dotted line: collisional loss.) 

For > a 0^2, the 



behavior. The slope of the latter decreases as mass grows 
slope is given by the BH formula (|4.24|) . 

Our main results for a heavy quark produced in the plasma are qualitatively 



represented in Fig. |T^. The transition between the quadratic and linear regimes 
occurs at a length scale which is ~ mm{L* , L* , E / {ij,M)) . For <^ agV ET^ it is 
L*, for asVET^ < < a.^^ it is L*, and for > a.E^ it is E/{^M). 

6. Concluding remarks 

The primary purpose of this review was, as for any other review, to bring together 
and present in a systematic way the known results scattered in the original papers. 
Another goal was to rederive those results using simple physical arguments rather 
than coming to grips with complicated formalisms. We tried, for instance, to explain 
in a semi-heuristic way the origin of the quadratic dependence AE'^ad ~ for thin 
plasmas or of the law Ai?rad ~ L\fE for thick plasmas. However, besides reviewing 
known results, we have made some new observations. 

First, we found that the law, which was generally believed to be a specific 
feature of QCD, is valid also in the Abelian case. The extra suppression at small 
L compared to a linear behavior ^E^a.d{L) ~ L is always present when a particle is 
created within the medium in a hard process. It is due to the fact that a newborn 
particle needs time to grow its radiation field coat and acquire the capacity to radiate. 
We stress that, although the Abelian and non- Abelian physical pictures and results 
are similar as far as the average radiative energy loss is concerned, the spectra of 
emitted photons in a QED plasma and of emitted gluons in a QGP are different 
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(for thick plasmas, the corresponding spectra are given in ( 3.30 ) and ( 4.20 )). The 
difference is due to different kinematics of photon and gluon bremsstrahlung, as 
discussed in sections ^ and ^. In QCD the presence of the extra graph of Fig. ^ 
(together with the graphs of Fig. ^ with appropriate color factors) broadens the 
gluon angular spectrum, see ( |4.6| ) and (^TTj). 

Another task we tried to accomplish was the systematic analysis of the radiative 
energy loss of a massive particle in the different regions of M and L. We have 
emphasized (see Fig. |15D that the mass effects play no role at small enough length, 
however large the quark mass is. For <^ asvET^-, there is no effect also at large L. 
For larger masses (one should distinguish the main regions asvET^ -C -C a.^E'^ 
and ^ OisE"^) the heavy quark loss A£'rad(-Z^) starts to deviate from the light 
quark loss at some critical length which decreases when M increases. 

Let us now discuss the question of the physical meaning of the mean free path A. 
The definition ( p.9|) relates A to the so-called anomalous damping ( of ultrarelativistic 
collective excitations with quark or gluon quantum numbers [35]. The anomalous 
damping depends on the total cross section {( ~ ncr*°*), rather than on the transport 
cross section ( |2.10| ). It is the latter rather than the former which determines the 



scale of different transport phenomena and a legitimate question is whether ( (or 
equivalently A) is a physical observable quantity. 

This question was studied in Ref. [50]. No way to measure ( was found there in 
the ultrarelativistic plasma, even in a thought experiment, but it was found that, in 
a nonrelativistic (Boltzmann) plasma and in a certain range of parameters, ( shows 
up in the argument of Coulomb logarithms describing transport phenomena. 

Coming back to the energy loss problem, we see that in most formulae A enters 
not alone, but in the transport coefficient 

q = ^ - al^T\ (6.1) 

For example, the estimate (|3.26|) for the electron radiative loss in a thick plasma 
is represented as A£'i.ad(-Z^ ^ L*) ~ aL^/qE. The parameter q describes how the 
transverse momentum of the particle grows with distance, {q\{L)) = qL. To be 
more precise, q is of the form (compare to 



/•|i|max rprp 

qr^n ^\t\d\t\oca'THn^, (6.2) 

J ^2 at jjL^ 

and thus depends logarithmically on the energy of the incoming particle. Recently, 
the coefficient of the logarithm in (|6.2| ) was evaluated analytically [51]. 

Seemingly, the parameter A may show up as a scale that distinguishes a very thin 
plasma L <^ A where the particle undergoes at most one scattering and the interme- 
diate region A ^ L ^ L*, where the multiple scattering kinematics is effective. In 
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our whole discussion, we indeed made a clear distinction between these two regions 
and treated them differently. 

As far as the radiative energy loss of a particle created in a plasma is concerned, 
nothing essential happens at L ~ A, as can be qualitatively seen on Fig. 12. However, 
in the case of a light quark produced in a QGP and at L <^ A, the parameter A 
enters the argument of the logarithm in ( |5.21| ) (use /i^ = qX). The region L <^ X 
is, however, rather marginal - it corresponds to a quark produced near the edge of 
the plasma. Since it seems unrealistic to "tag" such quarks, the observability of A 
in this case is questionable. In a more realistic case, A ^ L ^ L*, the logarithm in 
the large u ^ uJc contribution (|5.27| ) depends only on the combination ( |6.1|) . As was 
already mentioned, the small uj < Uc contribution may provide a logarithmic factor 
~ ln(L/A), and the scale ~ A can in this case be observed through weak logarithmic 
dependence. More studies of this delicate issue are required. 

For an asymptotic particle, the situation looks different. As is clear from Fig. ^ 
the dependence is essentially modified at L ~ A. In addition, the slope of the curve 
in the BH region L <^ A is given by the estimate ( |4.9| ) which explicitly involves A. 
On the other hand, it is difficult to imagine how a plasma (in thermal equilibrium) 
of size L ~ A or less could be created, as we already mentioned in the footnote at 
the beginning of section 3.1. 

Another attempt to pinpoint an explicit dependence on A is associated with the 
estimate ( 4.24| ) for the radiative energy loss of a heavy quark. We have seen that, 
when the mass is large enough, ^ o^gE'^, this estimate is valid not only for 
unphysically thin, L <^ A, but also for thick plasmas, see the dashed curve corre- 
sponding to ^ q^sE"^ in Fig- 0- The expression ( 4.24 ) involves the combination 
^/X = \/q/X. However, ( [4.24]) describes only the radiative energy loss. And as 
we have seen, when ^ cusE"^, the radiative loss is suppressed compared to the 
collisional one. For light quarks, the radiative and coUisional losses have different 
patterns: the characteristic energy of the radiated gluons is of order E, which is 
much larger than the characteristic energy transfer in one elastic collision. But the 
radiation spectrum of heavy quarks is soft. The spectrum is cut off beyond the scale 
fiE/M, which for large masses ^ dsE'^ is smaller than the plasma tempera- 
ture T. In other words, it seems impossible to separate for so heavy quarks the 
radiative component of the net drag force dE/dx and access /i/A and thus A. Quite 
curiously, for smaller masses, when radiative losses dominate, their value is not sen- 
sitive to A. For example, a nontrivial estimate ( |4.31 ) in the intermediate mass region 
asV ET^ <^ <^ OgE"^ depends only on q. 

Still the parameter A (and not only the combination ( |6.1|) ) seems to have an 
independent physical meaning. In fact, this parameter appears under the logarithm 
in the more refined estimates ( |3.29|) , ( |4.19|) for the light particle radiative losses 
(recall that T ~ A/i^ = qX"^). These estimates take into account the behavior ~ 
(A^ In A^) rather than just Nfi"^ for the characteristic effective scattering momentum 
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transfer. In the developed LPM regime, N = L*/X = ^/WJiXiF). Assuming that the 
estimates ( p.29| ), ([4.19|) are correct, the situation is then analogous to that observed 
in Ref. [50] for Boltzmann plasmas: the parameter A affects observable quantities in 
a weak logarithmic way. 

In other words, the physical status of A (or Q in ultrarelativistic plasmas remains 
unclear. But it is indisputably very useful as a theoretical instrument allowing one 
to obtain meaningful physical results for radiation spectra and energy losses. 

As a dessert, let us mention the fascinating issue of energy losses in A/" = 4 
supersymmetric Yang-Mills (SYM) theory. At weak coupling, there is not much 
difference with QCD, and we expect the estimates quoted in this paper to apply 
also to A/" = 4 SYM. The main interest of the A/" = 4 SYM theory is that in the 
large Nc limit, many quantities can be evaluated also at strong 't Hooft coupling^^ 
A = g'^Nc ^ 1, using the duality conjecture [52]. In particular, the drag force dE/dx 
acting on a heavy quark moving through a thick A/" = 4 SYM plasma reads [53] 



This estimate is valid when M ^ {XTE'^Y^^ [48,54]. The dependence (|6.3|) resembles 
the perturbative result ( [4.24| ). One difference is that the latter is valid in a different 
mass region, namely M 3> gE, and that it describes only the radiative energy loss 
which happens to be dominated in this region by the collisional loss. 

In Refs. [55] (see also [56]), the energy losses of light partons in a strongly coupled 
A/" = 4 plasma were estimated. The dependence 

— ~ -X^^E^tY' (6.4) 
dx 

for the mean drag force (for light partons, this quantity involves large fluctuations) 
was found. This is different from the perturbative dependence dE/dx oc V^- More 
studies in this direction are desirable. 
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They were obtained in a model where the particle is scattered on a set of static Coulomb 



sources separated by the average distance A. The presence of the factor cx Vln E in (3.29) is a 



robust model- independent feature, the origin of the logarithm being the same as in ( |6.2| ). On the 
other hand, in what particular way the argument of the logarithm is formed, whether it is E/T, 
ET/ii^ ^ E/{aT), or some other ratio, is a more delicate and difficult question. Only an exact 
model- independent calculation (possibly using the formalism of Ref. [37]) could resolve it. 
•^^Do not confuse it with the mean free path! 
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A. Typical momentum broadening in Coulomb rescattering 

Here we consider a charged (colored) particle with E oo moving in a perturbative 
QED (QCD) plasma, and undergoing n successive Coulomb scatterings. The range 

of the Coulomb potential is assumed to be much smaller than the mean free path 
A between two successive scatterings, so that the elastic Coulomb rescatterings are 
independent. We calculate the typical transverse momentum gtyp(^) of the particle 
after N scatterings. In the case of fixed coupling (QED), the result was found in 
Ref. [14] to be qly^{N) ~ fi'^NlnN. We present an alternative derivation of this 
result and generalize it to running coupling (QCD). 

Consider first the Abelian fixed coupling case. Coulomb scattering with trans- 
verse momentum exchange is associated to the normalized probability density 



The average momentum exchange {q"^) in a single Coulomb scattering is logarithmi- 
cally divergent. (The divergence is cut-off by the kinematical constraint on the maxi- 
mal transverse exchange |(?|max, but we focus on the E ^ oo limit where |g|max ^ oo.) 
On the other hand, the typical transverse momentum transfer gtyp, defined as the 
transfer such that the probability to have |g| < gtyp is 1/2, is well-defined. Solving 
the equation 



we easily find that gtyp in one scattering equals the Debye mass /x. 

We now determine the typical transfer gtyp(^) after scatterings, defined by 




(A.l) 




(A.2) 




2 



1 



(A.3) 



Representing the 5-function as 




(A.4) 



we obtain from (|A.3|) , 




(A.5) 



where 




(A.6) 
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The g-integral in ( |A.5|) can be done exactly, leading to 



1 roo POO a 

2=y^ drq,y^{N)J,{q,y^{N)r)[rKr{r)f = - rfr [riTi^f - Jo(gtyp(iV) r) . 

(A.7) 

In the latter equation and in the following, r is expressed in units of and gtyp(^) 
in units of /i. Integrating by parts and using {rKi{r))' = —rKQ{r) we arrive at 



N I drrJoiqtypiN) r) [rK,{r)f-' Ko{r) . 



(A.8) 



The equation ( |A.8| ) for gtyp(A^) has been derived from ( |A.3| ) without any appro- 
ximation. We now assume a large number of scatterings, ^ 1, and derive the 
asymptotic behavior of qtyp{N) in this limit. 

Clearly, when ^ 1, the integral in ( |A.8| ) is saturated by r <^ 1. We can thus 
approximate 



Knir) ~ — Inr 



r\ 1 
4 



rKi (r) ~ 1 In — ~ exp 



, 1 

— rl^^ 

4 r2 



Since A'o(^) is a slowly varying function for r <^ 1, we obtain from (|A.8| ), 



A^(lnl/r) / dr r Jo{qtyp{N) r) exp 



Nr^ , 1 
T- 1^ 

4 



The integral is dominated by the region 



A^ In 



1 



1 



Using this, we can rewrite (|A.1CI|) as 



A^lnAT 



l^A^lnA^ / dr r Jo{qtyp{N) r) exp 
Jo 



-{NlnN) 



4 



Introducing u = (NlnN) r^, this becomes 



du f gtyp(A^) 



2 VVA^lnA^ 



u e 



-u/4 



Using 



POO 

/ du Jo(Cv^)e-"/^ = 4e-^' , 
we finally obtain (reintroducing the dimension of gtyp(A^)) 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



(A. 14) 



qlJN)c^{\n2)-iN\nN)-fx\ 



(A.15) 
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This result immediately follows also from the expression 

1 ( <il 



for the probability distribution of the transverse momentum transfer after 
scatterings, derived previously in Ref. [14]. Indeed, defining the typical transfer as 
in ( |A.3| ), namely 

ciV/(g^A^)e(g,2,p(A^)-g2) = i, (A.17) 

and using ( |A.16D , we recover ( |A.15| ). 

The derivation above was performed for fixed coupling, and the estimate (|A.15 ) is 
thus valid for QED. In QCD, the running of the coupling should be taken into account 
and this brings about certain modifications. The effective coupling constant depends 
on the transverse momentum transfer . The normalized probability density of a 
single Coulomb scattering is now 

/i^ / ^ \ 1 ^2 In^T?^ 



Pf^l I -sVi I pi > ^ ^^QCD fA 18) 

n^JlQCD- + ^Q^J -^{q^ + ^^-Y In^^ ' ^""-'"^ 

where F{x) is a smooth function that tends to unity in the limit /i ^ Aqcd we are 
interested in, and which we have thus been allowed to set to F = 1. The analysis is 
done along the same lines as in QED. 

We first note that the expression ( |A.8| ) can be rewritten for a general scattering 
probability density P{r) as 



-N I drJoiqtypiN)r) Pir) . (A.19) 



N-l 



dP{r) 



When A^ ^ 1, we have typically r -C 1 (r being expressed in units of yU ^), implying 
P{r) ~ 1 (this can be easily checked a posteriori). Thus, ( |A.iy| ) can be approximated 
by 

\-N rfrJo(gtyp(Ar)r)exp[-(Ar-l)(l-P(rO)]|:(l-P(f)) . (A.20) 
Using ( [A. 18 ), we obtain (expressing q in units of /i) 



1 - P(r) = j £q P(g-)|QCD (1 - e-^"-^") ^ ^ 



lA^ 1 ^ 2 2 

(g2 + 1)2 ^ 4 



'QCD 



^^QCD 
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Comparing to ( |A.9|) , we see that the running of as displays itself in the factor 
as(/i^/r^)/as(/i^). Using now ( |A.11|) , we infer that the running of the coupling 
modifies the fixed coupling estimate ( [A.15| ) to 

It is interesting to mention that the typical momentum transfer at large N, 
in QED and QCD (see (|A.15| ) and (|A.22|) ), can be heuristically obtained from the 
following formulae 



<p(iV) - fi'N ^ (QED) , (A.23) 



lUN) f''''''^^ (QCD). (A.24) 

Finally, let us remark that ( A. 22 ) can also be represented as 



In In AT 



In-r^ + lnAT 

^^QCD 



Thus, for very large N , namely In ^ In ^^^^ , we obtain 



The scale is nothing but the average momentum transfer (g^) associated to the 
QCD probability density ( |A.18| ). In QED, this quantity is logarithmically divergent. 



but in QCD the running of as makes the integral for {q]) convergent even when the 
upper bound is put to infinity, as can be seen from ( |A.24| ). 

B. LPM effect and Feynman diagrams 

In the main body of the paper, we have operated mostly with heuristic arguments 
based on formation length estimates and single scattering formulae. The same re- 
sults can be derived by calculating the Feynman diagrams describing photon (gluon) 
radiation in the process of multiple scattering of a fast particle in the plasma. In 
this Appendix, we will not attempt to perform a complete diagrammatic analysis, 
but will present some illustrative calculations which might help to understand better 
the origin of LPM suppression. We will restrict ourselves to the Abelian case and 
mostly follow the discussion of Ref. [41]. 

We adopt, as we did in Appendix A, the model where a scalar massless particle 
is scattered on static centers with a Yukawa potential. 



^^^^^ np rf* . 
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Figure 16: Electron elastic scattering on two centers. 

where Xj is the position of the i-th center. Consider the case of only two such centers 
and assume Xi = 0, X2 = (x2±, z). Then the elastic scattering amplitude (see Fig. 0) 
reads 

Me, oc / ,^,f%fS\ + ^2 - g)e-^^-- . (B.2) 

The total momentum transfer q and intermediate electron four-momentum pi are 
given by 

q = P2-po ; Pi = {E,qi^,E + qi\i) , (B.3) 

where we have chosen po = {E,0±, E). In the model of static centers, the energy 
transfer in each elastic scattering stricly vanishes, Q'l = 5'° = 0, implying Pq = Pi = 

pI = e. 

We integrate over dqi\\dq2\\5[qi\\ + ^211 — q\\) by closing the contour in the upper 
gi||-plane and picking up the contribution of the pole of p\ at the value^^ 

PHI =E + q,\\^E-^. (B.4) 
Using the on-shell condition = we get q\\ ~ —q\/{2E) and we obtain from ( |B.2| ) 



where 



$scatt = z{pi\\ - P2\\) ~ ^ + g2±)^ - qi±\ ■ (B.6) 



■^^The contributions of the poles of + fj,^ and q2 + are suppressed by ~ e ^ 1. Indeed, 
the distance between successive scattering centers is z ~ A where the mean free path A satisfies 



A ^ l/(e-^T) ^ fi ^ ^ l/{eT) in a perturbative plasma. 
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Figure 17: Amplitude for plioton radiation induced by double elastic scattering 



To evaluate the elastic cross section, we fix the longitudinal distance z between 
the scattering centers and average |A1eiP over X2±. Integrating further over d^gl 
yields 

Consider now the process where the scattered fast electron emits an additional 
photon. There are three graphs depicted in Fig. 117. 



It is convenient to define the momenta Pi as in the case of elastic scattering, 
see for instance (p.3| ), so that Pi "do not know" about the emitted photon. On the 
other hand, the final momentum is now P2 — k rather than p2, and the intermediate 
momentum pi — k appears in the graphs of Fig. p!7| a and This shift of momenta 
brings a modification of the phase factors in the amplitude. For different graphs, 
this modification is different, and now one cannot suppress the phase factors as we 
did in the elastic case, when going from ( |B.5| ) to ( [B.7] ). 



Let us see how it works. Consider first the graph of Fig. [T^. The conditions 
(pi — fc)^ = and (p2 — = imply, respectively, 

Pi|l-i5-|t-^. Pn^E-'^-^, (B.8) 

where 61^2 are the angles between the direction of the emitted photon and pi (^2)- 
We assumed that the angles are small. Substituting the expressions ( |B.8|) into the 
phase $ = z{j)i\\ — P2\\), we obtain for the graph of Fig. p!7^, 

= $scatt + ^{01 - el) , (B.9) 



33 



Indeed, the radiation probability is dominated by small angles. 
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where $scatt is given in ( p.6|) . The corresponding contribution to the radiation am- 
phtude is 



-mS = exp i^^iel - el) \ A^e. , (B.IO) 



where = k±/uj and e'is the photon polarization vector. 



For the graph of Fig. 17c describing the emission from the final line, the structure 



-^0 ■ ^/^o transformed into ^2 ■ ^/^i- The phase factor is different from that in the 



graph of Fig. p!7| a due to different kinematics. Putting the intermediate momentum 
on mass shell gives in this case the condition Pi = rather than {pi — k)"^ = 0, so 
that the expression for is not modified compared to the elastic scattering case. 
We have 

$i?d = $scatt + f^2^ (B.ll) 

and 



M^i^'-^e.p\':^el}M,. (B.12) 



'2 



The graph of Fig. [T7| b provides two different contributions from the poles Pi = 
and {pi — ky = 0. They both involve the structure ^1 ■ e/ef, but the residues 
have opposite signs. In addition, the phase factors for these two contributions are 
different. For the pole pi = 0, the phase coincides with (|B.11D, whereas for the pole 



(pi — ky = 0, it coincides with (|B.9| ). The sum of all contributions can be expressed 
as A^rad = eA4ei£- J, where 

J = J,exp[z'^{el-el)] + J2exp[z'^el] (B.13) 



Ui Uq fc/2 ^1 



Each term in the sum ( p.l3|) corresponds to the radiation induced by elastic scat- 
tering on the associated center. The phase difference ujze\/2 between the two terms 
can be interpreted as the phase acquired by the photon of energy uj in the frame 
moving with the fast particle (see (|3.20| )). 

The result ( |B.13| ) can be easily generalized to the case of scattering centers. 
Let us assume that their longitudinal positions are Zn = {n — 1)A, n = 1,. . . ,N. 
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Then 



J{N) 



exp < t 



.u;X{N-l) 



. N 

^ n=l 



J = — - 



6'n 6n-l 



dn = Oq — 



771=1 



Af-1 



Y^ei (n=l,...,iV-l) ; $;v = 0, 



(B.15) 
(B.16) 
(B.17) 



The radiation energy spectrum is 

dl 



a 

duj TT^ 



N 



n=l 



(B.18) 



This should be averaged over qn± with the weight 



n 



2^2 



(B.19) 



The analysis of the expression ( [B.18| ) confirms the physical picture outlined in 



section ^ In particular, (i) the characteristic total scattering angle is 91^^^ ~ N / E"^ 



(ii) the contributions of different scattering centers in ( [B.18|) are coherent (so that we 
are dealing in this case with one effective scattering) when uNX9^^^ ~ ^^-^^tot ^ 1; 
i.e., when L <^ a/ XE"^/ (cj/i^). The scale ^^/XE^JJujJj?) coincides with ij^^{uii) defined 
in ( p.25| ). At uj ^ E, the latter coincides with the characteristic in-medium formation 
length L* {i.e., the coherence length whence one photon of energy ~ £^ is emitted). 

The results ( |B.15| ) and ( |B.18| ) have been derived for an asymptotic particle. In 
case the particle is created in the medium, the radiation amplitude can be obtained 
from ( |B.15[ ) by treating the position of the first scattering center as the creation point, 
and by suppressing the contribution of the graph analogous to Fig. |T^ describing 
the emission from the initial line. Suppressing the irrelevant common phase factor 
in front of the sum in (|B.15|) , changing N N + 1 and the numeration 1 — > 0, etc., 
we derive 



J, 



creation+A'^ scatterings 



N 



e^*o ^J-Jn e'' 



(B.20) 



n=l 



For = 1, we reproduce the result i f).2\j . 

The medium-induced radiation spectrum is 



UJ- 



dl 

I 

duj 



induced 



N 



d% \J\ 



^0 



(B.21) 
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where the square of the first term in (p. 201 ) (corresponding to the vacuum contribu- 
tion) has been subtracted, as we did in ( ^.3|) . We obtain 

2" 



induced 



N 



N 



9, 



N 



n=l 



(B.22) 



which is convenient to represent as^^ 
dl 



duj 



induced 




~ a 


f d% 1 


N J 





/ d'Oo <^ 2-^ ■ ^ J„ [cos(<l>o - $„) - 1] 

^ I n=l 



N 



+ Jn-Jm [C0S($„ - $„) - 1] 



(B.23) 



It is manifest that the induced spectrum vanishes when L = 0, since this imphes 
$„ = for all n. 

Consider first the contribution of the first term of ( |B.23| ) and concentrate on a 
particular term in the sum Yin- seen from ( |B.17] ), the phase difference 

$n ~ ^0 does not depend on 9n and hence on qn±. We can thus average over qn± 
before doing the integral. Averaging over azimuthal directions gives (see ( |5.4D ) 



dn~l — <in±/ E 9n-l \ _ 



e 



OLi ■ (B.24) 



Averaging further over g^^ with the weight (|5.6|), we obtain 



(Jn) 



(B.25) 



Thus, we get the contribution 



dl 

duo 



1st term 



n=l 



d^Oo 



01 oi_M-i+f^yE^) 



n-l 



m=0 



(B.26) 



where, in each term of the sum, the averaging over qm± with m ^ n should still be 
performed. 

Suppose L <^ L* (otherwise the physics is the same as for the asymptotic particle, 
and there is no point to analyze ([B.23| ) instead of (P.18| )). Then 1 / {ujL) > 1/{EL) ^ 
Nfi'^/E'^. Thus, bearing in mind that 6"^ ~ 1/{ujL) (to be verified a posteriori) and 



I Q{m^n)- 



/i, we derive 



m ^ n ^ 6 , 



— 9n ^ 

■m 



ujL E'^ 



(B.27) 



34 



We use Y.n Jn = - B^IBl and / d^Q^ {l/Ol - l/Ol) - 0. 
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We arrive at the estimate ( ^.71) of the spectrum, giving (see 



dl 



induced 

~ a^^, (B.28) 

1st term 



which contributes to the average loss as 

Ai^llstte™ ~ (B.29) 

The integral is dominated by the angles of order 1 / {ooL) , indeed. 

Let us show now that, at L <^ L*, the contribution of the second term in (|B.23D is 
suppressed compared to ( p.28|) and (|B.29|) . The sum involves ~ iV^ terms. Consider 
one of these terms, say, the term with m = 1, n = N . The phase difference 



2 

k=l 



does not depend on g„x and we can average J„ over it as before. In addition, choosing 
6i rather than 6^ as an integration variable, we can observe that $iAr does not depend 
on gi_L, and we can also average Ji over the latter. We obtain 



dl 

duj 



! d^P gl -4^-1(1 -cos 'I'ljv) 



In contrast to the integral ( |B.26| ), this integral is not dominated by large angles with 
$iAr ~ 1. Indeed, assuming 9k — Oi, we are led to the integral 

32 



/ ^ [1 - cos{u;L9')] , (B.31) 



which is saturated by 



N^<^9^ <^^, (B.32) 
ujL 

and exhibits a logarithmic behavior. Expanding cos {ujL9'^) (we are allowed to do it 
in view of ( |B.32| )), multiplying the contribution of one term in the sum of the second 
term of ( [B.23| ) by A^^, and neglecting the logarithmic factor which is irrelevant in 
the present discussion, we obtain the estimate 



dl 

doj 



induced 

~"(9?). (B.33) 

2nd term 



Integrating this over uo gives the average loss 

Ai?l2ndterm ~ " ^ ~ j " ^^'^^^ 

This is indeed suppressed by ~ oJc/ E ~ {L/L*Y <^ 1 compared to the contribution 
( P:^ ) of the first term of (^). 
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